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Introduction 

The  equation  vie   are  dealing  with  is  relevant  to  quasigeo- 
strophic  motion  of  the  atmosphere.  It  is 

The  mathematical  formulation  of  the  atmospheric  motion  and 
derivation  of  the  above  equation  using  geostrophic  approxima- 
tions can  be  found  in  reference  [1], 

Equation  (1)  is  called  the  geostrophic  conservation  equa- 
tion where  d/dt  denotes  the  following  operator: 

dt   9t   ox  dy 
and 

A  is  the  Laplacian  operator  sVc^x^+^^Ay^,  and  k^  is  a  constant 
determined  physically  from  Coriolis  parameter  f,  acceleration 

of  gravity  g,  and  mean  height  of  the  atmosphere  h  ,  i.e., 

2    2/  P 

k  =  f  /ghQ.  We  call  (^-k^)  the  Helmholtzian  operator.  Equa- 
tions (1)  and  (I.b)  show  that  \|/(x,yjt)  is  a  stream  function  and 
u  and  V  represent  velocity  components  in  the  x  and  y  directions. 

It  has  been  shown  in  [2]  and  [3]  that  it  is  reasonable  to 
solve  the  above  equation  when  \|;  is  specified  initially  over  a 
rectangular  domain  G,  i.e.,  -l<x<l,  -l<y<l  and  also  on 
the  boundary  of  G  for  all  time,  i.e.,  0  <  t  <  T,  and  vrhen  ^\j<, 
which  denotes  vorticity,  is  specified  on  the  boundary  as  a  func- 
tion of  time,  only  v;hen  fluid  is  entering  the  region,  but  not 
when  it  is  leaving  the  region. 

Ch.  B.  Sensenlg  [2]  proved  an  existence  and  uniqueness 
theorem  for  the  above  problem  when  the  domain  is  a  half  plane. 
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e.g.,  X  >  0.  Hence  we  shall  solve  the  problem  for  the  square 
numerically,  i.e.,  we  shall  show  that  a  particular  finite  dif- 
ference scheme  provides  an  approximate  solution  of  the  above 
problem  accurate  vjithin  a  reasonable  error  under  the  assumption 
of  the  existence  of  a  sufficiently  smooth  solution. 

The  general  procedure  of  this  paper  is  divided  into  two 
parts:  development  of  the  error  formula  for  Helmholtzian  (Part 
I  -  chapters  1-5)  and  establishment  of  estimates  for  the  error 
Itself  and  its  difference  quotients  (part  II  -  chapters  6-12). 
One  of  the  interesting  features  of  the  convergence  proof,  which 
is  necessitated  by  the  non-linearity  of  the  differential  equa- 
tion (1),  is  that  we  have  to  estimate  up  to  the  boundary  the 
bounds  for  the  first  difference  quotients  of  the  solution  of 
the  finite  difference  analogue  of  the  Poisson  equation,  in 
terms  of  a  bound  for  the  right  hand  side. 

VJe  shall  start  the  first  part  by  introducing  in  chapter  1 
a  particular  finite  difference  scheme.  By  assuming  that  the 
true  solution  of  (1)  has  continuous  bounded  derivatives  up  to 
order  k-   and  that  the  i|.th  derivatives  satisfy  a  Lipschitz  condi- 
tion, we  compute  the  truncation  error  (denoted  by  T(\I/))  and 
find  that  it  is  of  order  b,  i.e.,  T(\!/)  =  0(h). 

In  chapter  2  we  analyze  step  by  step  the  growth  of  the 
error  for  the  Helmholtzian  assuming  that  initial  error  is  zero 
and  that  the  difference  equation  which  approximates  the  diffe- 
rential equation  is  solved  exactly.   Under  the  restrictions 

that  /^t  <  A.V[wax(|A/^x^|  +  IA/Ayf)  J  ^^'^  '^   <  VNq  (where  | 
denotes  the  solution  of  the  difference  equation,  T  —  the  total 
time  after  n  time  steps,  i.e.,  T  =  n/\t,  and  N  is  a  constant 
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3 
depending  on  the  bounds  of  the  derivatives  of  the  true  solu- 
tion) we  find  the  following  formula  governing  the  error  of  the 
Helmholtzian,  i.e.,  IZLhe^^^"^^®^"^^  I  -  ^o^  "  °  ,  in  which  e  de- 
signates the  error  e  =  \Ir-^. 

In  chapter  3  we  consider  the  initial  error  of  the  Helmholt- 
zian  (denoted  by  E^°')  to  be  present  and  by  similar  analysis  as 
in  chapter  2  we  develop  an  error  formula  for  the  Helmholtzian 
and  we  shall  conclude  that  if  the  initial  error  is  of  order 
0(h)  and  if  the  restrictions  of  the  previous  chapter  are  satis- 
fied the  error  of  the  Helmholtzian  is  of  the  same  order  as  that 
in  chapter  2. 

In  chapter  I4.  by  introducing  a  concept  of  "time  layer"  wo 
show  that  it  is  possible  to  extend  the  total  time  (denoted  by 
T)  to  T  <  k/N   (where  k  denotes  the  number  of  "time  layers") 
and  by  a  similar  error  analysis  as  before  we  find  a  verified 
error  formula  for  the  Helmholtzian. 

Chapter  5  briefly  shows  that  the  above  established  conver- 
gence is  subiject  to  the  requirement  that  the  round-off  be  of 
order  O(h^) . 

In  chapter  6  we  obtain  from  the  bound  for  the  Helmholtzian 
by  means  of  maximum  principle  and  auxiliary  function  a  bo\ind 
for  the  error  e. 

In  chapters  7-10  v;e  make  the  necessary  preparations  for 
establishing  the  estimates  for  the  difference  quotients  of  the 
error.  We  simplify  and  reduce  the  original  problem  to  two  se- 
perate  problems:   the  non-homogeneous  and  the  homogeneous  one. 
Solving  the  non-homogeneous  problem  vje  come  across  the  solution 
of  the  so-called  "discrete  potential  equation"  whose  properties 
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are  discussed  in  chapter  2  using  von  der  Pol's  and  A,  Stohr»s 
results  in  [5]  and  [6],  respectively. 

In  chapters  8  and  9  we  obtain  subsidiary  bounds  for  the 
solution  and  its  difference  quotients  of  the  non-homogeneous 
problem  with  the  help  of  the  boimds  for  the  solution  of  the 
"discrete  potential  equation"  given  by  A,  St'^hr  in  [6]. 

In  chapter  10  we  obtain  bounds  for  the  difference  quotients 
of  the  solution  of  the  homogeneous  problem  applying  the  method 
given  in  Tamarkin-Peller  [I|.j. 

Chapter  11  establishes  the  necessary  estimates  for  the  dif- 
ference quotients  of  the  error. 

In  chapter  12  vje  show  that  using  the  reflection  principle 
we  can  extend  the  above  mentioned  estimates  up  to  the  boundary. 
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PART  I 


1.  Truncation  Error 

Let  us  define  a  rectangular  grid  of  points  by  coordinates 

3^k  "^  ^^^  ^  ^   0,1,.. .,p 

(1.1)  ,. 

yji^  =  5Z\y  -^=  0,1, ...,q 

and  denote  quantities  at  the  point  (k,j?)  by  subscrpt  k, i^.   In 
the  same  manner  we  replace  t  by  m/\t .  m  having  integral  values. 
Let  us  choose  the  mesh  width  to  be  /\x  =  /'.y  =  h.  Prom  the 
following  diagram 


y  I 


k,/+l 

k-1,/  k*,/  k+i,/ 

k^ir-i 


we  approximate  the  derivatives  in  (1)  using  centered  differen- 
ces as  follows: 

,  .^  ^(x+A,x,y;t)-t(x-Ax,y;t)  /\     \+l,r'''k-l,  S^ 

^^  2^x  £^x^        2h 

^  ''^(^>y^Ziyi^^-''^^^>y-Ay;t)  _  A_^  _  '''k,g+i"^k,i?-i 
""  ■  ^  "  2Ay  "  Ay   '         2h 


(1.2) 


and  we  replace  the  Laplacian  of  quantity  G  by  the  following  fi- 
nite difference  approximation  /V.  vjhere 

hr 
The  differential  operator  d/dt  is  approximated  in  the  fol- 
lowing manner  taking  into  consideration  the  signs  of  u  and  v. 
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Case  1  If  u  >  0,  V  >  0  then 

^WT'^.^'X     ay's  ^^ 

il.k)  ^   u(x.y,t)[g^^>^>^^t^'"^'^^^'^3^ 


+  v(x,y,t)[ j^^ •^=i — ]  , 

Case  2   If  u  >  0,  v  <  0  then 

(1,5)  +  u(x,y,t)[ ]  + 


,  vu,y,t)[g"'-'-^"Ay.t)-6(^,y.t))  , 

h 
Case  3  If  u  <  0  and  v  >  0 

r  5   ,„■£>  ,,,  5,1-  g(x,y,t-f^t)-g(x,y,t) 

(1.6)  .  u(x,y.t)[^^^^^^-^f-g^->^-^^]. 

.v(x,y,t)[S(iLa.I):i?(iiIi^.ZiiIj  . 

Case  U  If  u  <  0  and  v  <  0 

!^u^+vj|lg  =  g(x,y,t+/\t)-g(x,y,t)  + 
(1.7)  .u(x.y.t)[«'^"^-'y-g'-'^'^'i. 


[  .  -  •• . .-.  I ., .,  ,_  1,/  ^  ,. , ..'   J- 
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Using  the  finite  difference  approximations  (1.2)  and  (1.3) 
equation  (1)  becomes: 

(1.8)  [^-(^^l')^+(^^)^](Ab^-k^^)  =  T!(^) 

At   Ay   A^   Ax   Ay      ^ 

where  /\//\^   denotes  the  finite  difference  approximation  to  t?/cbc, 

and  similarly 

3_  ::  A_     5  -^_ 

^y  ~  Ay    '  ^*  '  At 

and  T(il/)  denotes  the  truncation  error  (\|r  satisfies  the  finite 
difference  equation  except  for  the  truncation  error).  Prom  (1.6) 
we  shall  calculate  the  truncation  error.  According  to  (1.3) 

(1.9)  Ah^-k^*  =  ;^t^k+l,i^-^Vl,?-'*k,M-''^k,^-l-^^k^^-k^\j^ 

and  if  we  assume  that  the  fxonction  t(x,y,t)  has  bounded  deriva- 
tives up  to  order  i^.,  we  can  apply  Taylor's  expansion  and  get: 


A  ,,2,   If,  .A  ^"'id^  .  (A^)^  S''>M  ,  (A^)^y\it 
Ah*-''  ♦  =  ;^t*^,W  ^^  *  -1—  ^p-  *  -6—  ^;r  * 

c)X 
2    9x2      6    ^y3      2ii.     ^jK 

^  ,      A   ^*kj?  ,  (Ay)^  ^^^kfi     (Ay) 3  3\i>  ^ 
-  *,  r  Ay  —  - -T- -^  - -T- T;r  ^ 

In  our  case  Ax  =  Ay  =  h,  hence 


'r-c:^ 
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^   K   J?+l  !?  -^^k  P-1   Q.     2 

and  we  get 

(1.10)  Al,^-k2*  «  H(^)+  |j^k^(^) 
where 

(1.11)  H(^)  =  /\^-k^^ 

^^'^  ).-         L-         l±- 

(1.12)  Rkn(^)  = — ^^+ iLiSziiJL., kV-fl,^^   k,V  1,>^ 

^^       ax^        ax^      3y^       ^y^^ 

where  ^^?i,  t,+t  »/^x^  denotes  the  value  of  ^^\J//3x^  at  the  raid- 
point  of  the  segment  joining  the  points  (k,^)  and  (k+l,il),  si- 
milarly ^^?j.  j._T  g/<5^  ^h®  value  of  ^%/^x^   at  a  point  between 
(kj)  and  (k-l,i^),  as  well  as  ^^^k  J+i^)?/-^y^  and  ^^k,  ^-i  j/^y^ 
denote  the  values  of  ^H/c^  at  points  between  (k,P)  and  (k,i^+l) 
and  (k,J?-l)  and  (k,V),  respectively. 

Coming  back  to  (1.8)  and  substituting  (1.10),  we  get 

T(^)  =  [4r-(7^^)x"  +  (^^)^HH(lI/)+|r;Ri^„(^) 

A*  Liy    Zi^  Za^  /^y     2ii  Kx 

,1.13)  =  ^H(*).(^«^H(*)^(^t)^H(*). 

Let  us   first   evaluate    {/\//\t)E{i<) ,   i.e., 

At  At 

where  superscripts  m+1  and  m  denote  the  value  of  H  at  the  time 
steps  ra+1  and  m  respectively.  Having  in  mind  the  above  made  as- 
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sumption  we  see  that  E{^)    as  defined  by  (1,11)  has  bounded  de- 
rivatives up  to  order  2,  hence,  applying  Taylor's  expansion,  we 
shall  get; 

A.H(t)  =  -^[H^(vl')+AtH?(^)+  ■^^i?t(t)-H"(^'/)]  . 

At  At  *  2      ^t 

Considering  the  first  time  step,  i.e.,  letting  m  =  1 

(1.15)  ^H(^)  =  H^{M/)+  =|^Ht^(il/)  . 

Next  we  evaluate  A/A^^^^!')  ^^^   A/AyH(i}f)  expanding  it 
into  Taylor's  series  with  respect  t;o  x  and  y  respectively,  ob- 
taining 

A  Ax 

(1.16)  ^=^H(^)  =  H  (^)-  =4-^,(v|/) 

and 

(1.17)  ^H(^)  =  Hy(t)-  ^Hyy(^)  . 

However,  to  estimate  T(\!;)  vie  shall  need  some  more  estimates,  as 
(A/A^)^  and  (A/Ay)**   Consider 

^^  =  K^^l'K^'l  =   ^  (.,y^t)  =  .  (x,y,t).e  . 

A7  2Ay  ^  ^ 

Assuming  that  ^i     satisfies  a  Lipschitz  condition,  I.e., 

|\|/y(x,y"\t)-ty(x,y,t)|  <  k3_  I y'''"-y  I  <  k^^A^  • 
We  see  that  e  must  be  proportional  to  k^^ Ay ,    i.e., 

1^1  <  ^lAy   • 

Therefore 


(1.18)  TI^  ""  *y+^iAy 


A 


-9'. 
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and  similarly 
(1.19) 

Next  we  find 

A  rh' 


A 
Ax 


*  ~  ^v-'^^pA^  • 


Arh' 


A  rh^ 


T^t|-Ri^n(t)j  =  |r--^[R(x,y,t+At)-R(x,y,t)]  . 

At  2^  ^     ^^  At 

If  we  assume  that  Ri^j(^),  defined  by  (1,12),  satisfies  a 
Lipschitz  condition  v;ith  respect  to  all  the  variables  (i.e.,  the 
fourth  derivatives  of  ij/  satisfy  a  Lipschitz  condition) 

|R(x,y,t+/\t)-R(x,y,t)|  <  jj^^/\t 
|R(x+A,x,y,t)-R(x,y,t)|  <  H^A^ 

|R(x,y+A7»t)-R(^»y,t)|  <  i^Ay  • 

Considering  the  above  we  shall  get  the  following  estimates: 


(1.20) 
(1.21) 
(1.22) 


9^[|-R(^)] 

At  ^^- 

^[|rR(t)] 

Ax  2^ 

/'^   h2 
=^[|r-R(^)J 

Ay  2^ 


<  h^i?'^ 


2r"- 


<  h 


<  h^i'* 


3  • 


Denoting  the  following  parts  in  (1.13)  t>y  B  and  G 
(1.23) 


A      /   A  A   A 

B  =  ^=±-H(v[r)-(^^\|/)^=^H(vI;)+(^^\l/)=^H(^) 

At         Ay   Ax         Ax   Ay 


Ax   Ay  ^^ 


61 


X- 


Jl.f— -r^v. 


>.   -  \ 


-ii- 


j\ 


..■•»•.- 


bnil    aW  cTxol'i 


A 


;  ■. .  ^       .   a  ■■■  I  > 


rSli        .• 


;■■,■■■     -•      ■•!■■• ,  t 


SJlO-5  ^liY!i:ctfc^;      .^; 


;/.'  \  .: 


(  .'■**-'.■.. 


'U>'  o^^  ^l;.^ •.■/'/.. ^!.un»w 


I  \ 


U^.i; 


sdi   r<Mi:ffc)fi^;J 


■!■     i    ( 


•  I  ■■-•  .1  V  ; 
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from  (1.13)  we  have 

(1.25)  T(\!r)  *  B+C 

and  returning  to  the  original  expression  for  E{^) ,    as  indicated 

in  (1.11),  and  substituting  (1.15),  (1.16)  and  (1.17),  also 

(1.18)  and  (1.19),  into  (1.23),  remembering  that  /\x   s  ^y  =  h, 

we  get 


At 


B  =  (^^-k^\i/)^+  -^{/;:j-k^^)  ^^••{^^+^-^h)[i/\j/~i&^) 


2. . 

X 


-  |(Zi^-k^^)xx^-'(*x'"V^^^A!'-^^^)y-  |(A^-k2^)   ] 


At, 


B  =  (A^-i^^^)t-K(A^-^^^')x'^^x^A^-i^^^)^+  %-(A^-k2i!/)^ 


(1.26)    -  k^h(/N^ilr-k2,l,)^+i|/^(^^-k2v|;)  +  |l'y(A^-l^^) 


S  '-^  "XX 


k^h^  rr-     ^       o-      V.2 


+  -V-(A^-k^^)xx-  |^x(At-k^^)yy-k2|-(A^-k^^)yy  . 
As  \|; .  Is  bounded  and,  by  previous  assumption,  has  bounded 
derivatives  up  to  the  l+th  order,  vie   can  introduce  the  following 
bounds  for  certain  expressions  in  (1,26): 

(1.27)   k^l  -  ^1  J  I'^y'  -  ^2  '  l(A^-k^^)xl  ^  ^2   * 

|(/v,^-k2yl,)y|  <  N^ 


(1.28)   |(A^-l^^)ttl  ^^5  5  l(A^-k'^*)xxl  -^^6  5 


|(/^^^.k'^ijr)^y|  <  N^  . 

Considering  (1)  and  (l.a)  and  the  bounds  in  (1.27)  and 
(1.28),  we  get  the  following  estimate  for  (1.26): 

At  V,      ^il'i^    V,       1,2 

|b|  <  s^N^+k^hN2+k2hN.  +  ^'^2^^+   -^N^+  |N^N^+k2%-N^ 

or 

(1,29)  |bI  <  M^At+M^h+I^h^ 


4    :'Ui.; 


,vIX.JC)    .= 


■  p.a»'>  ^. 


T  p.  i  f  ■ 


:-r^ 


.    ...(^!-"-v 


-'-si/:.^? 


\U^      M'''      ?V'->>Av    ■' 


•I  -f-^t  •*-i 


i  _  J 


V'- 


,.ur^>i-!-^,)..lj... 
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where 
(1.30)  M^=-#;  M^  =  W^^aVi^*  2   •  ^2  =  ^  2  ^ 

Similarly,  using  (1.20),  (1.21)  and  (1.22),  (l.l8)  and  (1.19) 
and  bounds  for  \|/  and  ^   ,  \^e   get  the  following  estimate  for  C! 

/Xt  2i(.  1^'^       ^y    ^  2i^  k*.      £^   ^y 

•    t^^k^^^^l 
|C  I    <  £5[h2+h!/y+k^hl  f^h2+|4f^+k2h|  v|h2 

|C|    <  h^[.r£+N2y3+k3_htt+%^2'^IC2'2^^ 

(1.31) 

|G|    <  h2[r£+N2J^3+N-Li/|j+h3[k3_y'^'+k24'] 

|C|    <  Mgh^+M-h^    , 

And  finally  by  (1.25),  (1.29)  and  (1.31),  we  get  the  following 
bound  for  the  truncation  error  T(\|/) 

v,^  xmIv.'- XM  V.3 


I T  ( \|/ )  I  <  M^^t  +M^h+M2h'^  +}''l^hr  +M-h- 


or 


2, 


(1.32)  |T(i|/)|  <  M^At+%h+0(h'')  . 


2 •  Estimate  for  the  Helmlioltzian  (Initial  Error  E^°     =  0) 

VJe  shall  now  proceed  to  estimate  the  error  of  the  Helm- 
holtzian  after  the  first  time  step  assuming  that  the  initial 
error  is  zero.  We  will  designate  by  J  the  values  of  the  solu- 


1         •      ^f 


•  c... 


?!«.;.!■ 


•••  U  ^ 


n'i--^ti. 


-a 


l3n« 


V.r^-  .  J 


,      /'-.» 


'.   -   >     ^    ! 


I  r        M. 


*■   .      4.  J.   t  •.  — 

-III r-j xr»M    o. 
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tlon  of  the  difference  equation  which  replaces  (1)  and,  for  the 
time  being,  ^^re  shall  ignore  the  roimd-off  error  and  assume  that 
all  the  arithmetic  steps  involved  in  solving  the  difference 
equation  are  carried  out  with  infinite  precision. 
Hence 

(2.1)  [4_.(^|)4_+(^|)^](  A  |.k2|)  =  0  . 

A*   Ay   A^   A^   Ay      ^ 

Prom  (1.8),  we  have 

(2.2)  [4^_(^^)^+(4„^)^](A  ^r.k^^,)  =  T{^)    . 

At   Ay   A^    Ax   Ay 

If  we  add  to  both  sides  of  (2.2)  the  following  auxiliary  expres* 
sionj 

Ay   A^    A^    Ay      ^ 

and  rearrange  the  terms,  we  shall  get: 

f~-(~zI)~+(^I)~^HAh^-k^^)  =  T(i|;)+[(^»1.-^|) 

At    Ay   A^    Ax   Ay      ^  Ay   Ay 

A    A     \    / 
.  4±-+(^|.  ^^)4i-](A  ^.k2^)  . 

Ax   Ax     Ax    Ay 

By  subtracting  (2.1)  from  (2.3)  and  denoting  the  error  by 
e,  i.e.,  \^-^  =  e,  we  get 

f^-(f-D^-(f-I)^](A^a.lc2e)  =  T(^).[(^e)^  - 

At   Ay   Ax   Ax   Ay  Ay   Ax 

Ax    Ay      ^ 

or 


3iq  sonnilr  uo  t-siTtso   yT.e 


\ 


r-  X  ■''•  '•         /"^ 


V-j:  "i;^  ^-^i-  ■'•--'•  ^-^  '''J^ 


A. 


'\65   an';        '■  iU0'X4    i.'".':?    ^?i?-i: j"  -. '■/j:3    ^cI 


■  .^  -  n-v 


1  .u       ■"Z-i        V^ . 


\'- 


^0 


11+ 

,^  ^,    At     ""  Ay   A^   A^   Ay      ^ 

(2,5) 

Ay   Ax   A^   Ay 

and 

Z^t  Ay     A^   Ax     Ay 

(2.6)  (A,^^.k2vi,"i)+[(4-yi)A^.(4-if)4.] 

Ay     Ax   Ax     Ay 

(Ah^-k^e)^ 

where  the  superscripts  m+1  and  ra  denote,  as  before,  the  m+1 
and  the  m   time  steps,  respectively. 

For  the  first  time  step  we  have  to  take  m  =  0,  and  by  as- 
suming that  there  Is  no  initial  error,  i.e.,  e*^  =  0,  and 
(Av,e-k  ©)°  =  0,  as  well  as  =^e°  =  0;  ^^e®  =  0,  (2.6)  becomes: 

^  Ay  Ax 

(Ah^-k^e)"^  =  AtT(^) 
and 

l(Ah^-k^e)^l  <  At|T(t)l  . 

Using  the  estimate  for  T(\!/),  i.e.,  (1.32)  of  the  previous  chap- 
ter and  for  convenience  denoting  ( A.j-O-k  e)-*-  by  A,  e^  '-k  e^^ 
we  shall  get 

I ^e  ^  ^ ^  -k^e  ^  ^ ^  I  <  A* [ "o At+^%h+0 ( h^ )  ] 
(2.7) 

lA^e^'^-k^e^lM  <5(1) 
where 
(2.8)  6^^^  =  MQ(/\t)2+M^/\th+0(h3)  . 


^M^ 


\ 


■^■i^n 


.;.     ■ '^'^(:A^-.^,^,' 


bt:a 


•'Z-^r 


■+•■  ■■        )■ 


.*    'f!?;. 


s  fj  Smt:- 


tr: 


\  \ 


^'c.rcio 


(  r\ 


S:Zj'\fii   uW 


x ; 


rt 
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In  part  II  using  (2.7)  we  establish  the  following  esti- 
mates for  e^-'-^  and  its  difference  quotients  ^-e^-*-^,  =-e^-'-^: 

A^         Ay 

(2.9)  le^^M  <  6^^^  ,   I^^^M  <C  e^^Mloghl  , 

Ax      "  ° 

l^e^^M  <  C  e^l^log  h| 

Ay  "" 

which  we  shall  presently  use  in  the  follovjing  error  analysis. 
Let  us  now  consider  the  error  of  the  Helmholtzian  after 
the  second  time  step,  i.e.,  xn  =  1.  Prom  (2.6)  we  shall  then 
have  the  following  expression: 

A.e(2).k2e(2)  =  ^tT(.)+At[  (^e^^)  )^-(4-e(l)  )4-] 

Ay       A^   Ax       Ay 

(2.10)  (Aj;(l).k2^,(l))H-[i+/^t(^l(l))^. 

Ay       Ax 

./Xt(^l(l))4.j(/^e(l)-k2e(l))  . 

Ax        Ay 

In  (2.10)  let  us  denote  by  K^"'"^  and  L^-""^  the  following 
parts: 

k(1)    =  ^t(4_e(l))^(Aj,^(l)-k2^(l))./:,t(4_e(l)) 

Ay       Ax  Ax 

(2.11) 


Ay     ^ 

(2.12)    L^^)  =  [i+/xti^t^h^-/:\,t{^-^^h^~]{AJ^^--'^^^^h 

Ay        Ax        Z^        Ay 

(2.13)  Ah®^^^-^^®^^^    =  AtT(\!')+K^^^+L^^^ 

and 

(2.li+)         lAh^^^^-^^^^^^l   ^  At|T(t)l  +  |K^^M  +  |L(lM    . 


■C'.i 


~.r:ip,'^     -■,—•;  ()-rr->[' .1' n*?'     pf-' 


rj  TI 


jn    .-■•i.;,. 


U^.^ 


^^S^ 


\ :.'  ^,,r:!  (i)?»   r.  :,   ^i 


l" 


v;..-:. 


fCfiSi 


a^ 


f(i^)^.i% 


'  .'.  •  "* 


U')r 


•T  r- 


**     '  -^  '■ ,  r  '^     • '  "f 


(X) 


}T 


^cA 


ijjq 


f  f.  <• .  c  ■, 


"r!.ZJi>! 
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Let  us  first  get  the  upper  bound  for  L^  '  and  later  for 
K^-*-^,  Prom  (2.12)  we  have 

Ay   A^ 
Ax        Ay 

and  if 

_  A_-r(l)  =  u  >  0 

Ay 

(2.15) 

^=^1^^^  =  V  >  0 

we  take  the  difference  quotients  as  proposed  in  (l.ij.)  and  ob- 
tain: 

(2.16)      -  (A,a'l)-.2e(l'),.,^,)-4_i(l))4l[,y^^e'^) 

^   Ay  ^   Ax 

^  Ay  '     ^   A^ 

(/\h®   -^^  ®   ^k,Jl-l  • 
Next  we  shall  shox^  that  =^1)^  '  and  -=-TI^  '  are  bounded. 

Ay  Ax 

Recalling  that 

|(1)  =  ^(l).e(l) 

we  get 


'•^A-i 


X 


Hi^ 


'a-'jl-' -'s^,_ 


(\r. 


>  :-U 
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(2.18)  ^^  ^^  ^^' 

,4_|(1)|  <,4_,(i),.|4.e(i)| 

Ax  A^  Ay 

and  similarly 

(2.19)  |4_|(1)|    <   |A.,^(l)|,.|A.e^lM    . 

Ay  Ay  Ay 

Prom  (1.18)    and   (1.19)   we  obtain: 

(2.20)  |4_^(1)|    <    Ki/^^M+kph 

A^  -       X  2 

(2.21)  l^t^^M    <    ki^M+k,h   . 

Ay  ^ 

A  A 

At   this  point  we   shall  not  use   the  bounds    for  —^e  and  =^0 

A^         Ay 

mentioned  in  the  beginning  of  this  chapter  in  (2.9)  but  evaluate 
==i-e  and  ===-6  as  follows: 

Ax         Ay 

A  e(l)      -e(l^     . 

A  p(l)    _  ^Ic+l.r^k-l,^. 

Ax       "         2^ 

1^7(^)1    <^(|e(i^I.IeUl,,|). 
Using  the  bound  for  e^    '    from   (2.9)   we  obtain: 

'^^  1-211  h 

and  using   (2.8)   we  get 

(2.22)  l^e^^M    <  0(h)    . 

A:^ 

Similarly,  we  would  get 

(2.23)  l^^^^l    <  0(h)    . 

Ay 


tj: 


(n,../^j.ru.A-.  ('-'- 


^a' 


f  r 


■{,i'-isL'    .        --fiO 


'.V;.r  {eX.I}   t'n-.-   (8x.i; 


o» 


S) 


aiiU 

:fA 

ri;t  nJL   ft 

■    .       ■;-( 

.   „-4   ; 

xA 

/\, 


•   ::'JU 


s^  n 


'i.' 


(es.sj 
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Combining  (2.18),  (2.20),  (2.22)  and  the  bound  for  \!'  from 
(1.27) t   we  obtain: 
(2.21;)  l^-I^^M  <  NT+0(h)  . 

Similarly,  combining  (2.19),  (2.21),  (2.23)  and  the  bound 
for  \1/  from  (1.27),  we  get: 

(2.2^)  |^|(1)|    <  Np+0(h)    . 

Ay 

At 

If  in   (2.17)  we  denote  ~~  by  X  and  impose  the  following 
restriction 

(2.26)  l-X[(-^l^^b-(^'^l^^h]    >0 

Ax  Ay 

or 

(2.27)  A  < 


max[|^l^l)|.|#-I^'^n 

Ax  Ay 


then,  together  with  (2.15),  (2.21+)  and  (2.25),  the  above  rest- 
riction implies  that  (2.17)  is  an  average  formed  with  non-nega- 
tive weights  which  add  up  to  one.  Therefore, 

and  consequently,  since  by  (2.7)  I  Ah®  ^■''^-^^®^'''^  1  -  ^^"""^  » 
(2.28)  IL^I^I  <  5(1)  . 

If  (2.15)  is  not  satisfied,  which  means  that  coefficients 
of  (2.17)  are  of  different  signs,  we  can  always  adjust  the  dif- 
ference quotients  talcing  either  forward  or  backward  differences, 
whichever  is  necessary  to  get  a  convex  combination  for  (2.17) • 
In  all  cases  the  restriction  for  A  will  be  the  sam.e,  which  is 

(2.27).  If  we  denote  Ah®^"*"^"^^®  ^"""^  ^y  ^kV »  *^®  above  can  be 
contracted  as  follows: 


•>        .()     >r   rt 


^r:uo.•f  f.-^:t   f>;v^    '^-c,  -^ 


V       - »    .  / 


Sni 


rv^r 


il 

(Ob.. 


~<-,itiii  = 


.i. 


i"S 


yu 


<-C) 


o   ?= 


OCI^ 


.^tb   er::!    ir^r'l 


^.3  0pnc'- 


t    \    '.      -  i 


a ,?:    r-,  r>  ■ 


Ay  Zix  ' 

.^ia)=u<o    Implies    ^iV=Milr=;il'^ 

Ax  Ay 

^i'l)  =  V  <  0   Implies   4^^V  =  H]ln-°i'^^   ■ 

Ax  Ay 

Then  (2.17)  becomes: 

(2.29)  l(1^  =  [l-A(|u|  +  |v|)^V+^iu|Gk_3i5n  U, J^-^^'^'^i^-sign  V 
where 

[  -1  if  g  <  0 
sign  g  =  I     0  If  g  =  0   . 
I  +1  if  g  >  0 

To  get  an  upper  bound  for  K^  '  we  shall  use  the  estimates 

for  e^-*-^,  =^^-^^,  — e^-'-^  vjhich  will  be  established  in  part  II 

Ax         Ay 

and  which  were  briefly  stated  in  (2.9).  They  are  of  sufficient 
order  to  insure  the  convergence  of  the  proposed  finite  differ- 
ence scheme.  By  (2.11) 

(2.30)  IK^^M  <  /s^t\^l-e^^)\\^i/K^y^Ki,^^^^))\^ 

Ay         Ax 


Ax 

Prom  (1.10)  we  have 


Ax         Ay 


2...  _   „.,...h2 


Ah^-k  >!'  =  H(\!/)+|^Ri,5(^) 

where  H(\l/)  and  \d^)   were  defined  by  (1.11)  and  (1,12),  res- 
pectively. 


J) 


n?4 


fi 


i     ! 


eT-Jl^ 


.*  V ' 
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Hence 
(2.31)  ^-(Ah'J'-^^^^)   =  ~-H(^)+^t|rW^^^ 

Ax     ^  Ax        Ax  2^  ^ 

and 

(2.32)         l^( Ah'J'-^^^^')  i  <  l~-H(t)l  +  l^[^Rka(^^^ 
Ax  Ax  Ax  ^-^ 

but  by    (1.11)    and    (1.16) 

^H(^!/)    =    (A*-k^^)x-|(AH'-k^t)xx 

Ax 

and  using  corresponding  bounds  from  (1.27)  and  (1.28)  we  obtain: 

(2.33)  |™-H(v[r)|  <  N  +|n^  . 

Ax 

Substituting  (2.33)  and  (1.21)  into  (2,32)  we  obtain: 

Ax  :s   o  ^ 

(2.31^)  A 

I— (Ah^-k^^)i  ^  N  +0(h)  . 

Ax     "^  ^ 

Similarly  using  (1.1?),  (1.11)  and  (1.22)  and  the  corres- 
ponding bounds  from  (1.2?)  we  get: 

A:y 
(2.35) 

I— (Ah^-i^^^)  1  <  V°(^^  • 
Ay 

If  we  substitute  (2,3l4.),  (2.35)  and  the  necessary  bounds 
from  (2.9)  into  (2.30)  and  keep  in  mind  (2.8),  vie  shall  obtain: 

|K^-^M  <  A^C^e^^Mlog  h|[N^+0(h)]+AtCQ5(^^ 

•  I  log  h|[N^+0(h)] 


t  / 


\  ■■4-  t  -'■ 


'    -'"'tx^^    ■      Si 


irdBi 


,\c     ■.-■\ 


A    ^  / 
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(2.36)"        IK^^^  I  <  N^Ate^-^Mlog  h|+0(h3) 

where  N^  =  Cq(N-+N^), 

If  we  consider  (2.28),  (2.36)  and  reball  that  A*|T(\lr)|  < 
<  5^  '  then  by  (2.11|)  for  the  second  time  step  we  get: 

(2.37) 

lAh^^^^-k^e^^)]  <  25(^UAtS^^^N^|l0S  ^I  =  5^2) 

(2.38)  lAh^^^^-l^^e^^M  <  ^^^^  • 

It  Is  obvious  that  the  method  used  in  part  II  to  get  the 
previously  stated  bounds  for  e^"*"',  ==--e^^  and e^  '    from 

Ax     „        Ay 

(2.7),  can  be  applied  to  obtain  from  (2.3°)  the  estimates  for 
q(2.)  ^   A_q(2)  ^^^  =^q(2)^  That  is,  we  shall  obtain: 

Ax  Ay 

le(2)l  <5(2)  ,   |4.e(2)|  <c  5(2)|iogh( 

Ax 
(2.39) 

and   |^e(2)|  <  c  6(2)|iog  h|  . 

Ay  ° 

In  a  similar  way  as  before  from  (2.6)  we  shall  obtain: 
(2.I]-0)     lAy^^'^^^^^U    <  At|T(^!/)|  +  |K(2)|^|L(2)| 

where  this  time  we  shall  obtain  the  following  corresponding 
bounds  for  K^^^  and  L^^^ 

(2.1+1)         IK^^^I  <  W^Ate^^Mlog  h|+0(h^) 

and 

(2.1^2)         |L^2),  <,,ax|Ahe^^^-k^e^^^l  <S^^^  • 

If  we   combine    (2.I|-0),    (2.!.|.1),    (2.U2)    and  recall  that 
At|T(\l/)|    <  e^"*"'   we  obtain 


TO" 


ECO   j.:.r.. 


S) 


^H»S) 
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If  we  substitute  Into  (2.[|3)  for  5^  '  its  value  in  terms 
of  6^  ',  i.e.,  (2,37)  and  leave  out  the  superscript  for  6^"^' , 
i.e.,  e^"^'  =  5,  v;e  find  that 

lAh^^^^-k^e^^M  <  6+WQAt(25+^t5N^|log  hDIlog  h| 

+  25+/^t6N^ I  log  hi 

=  5(3) 

IAhO^^^-k^e(3)l  <5(3)  . 

Prom  (2.UIj.)  in  a  similar  way  as  before  v/e  shall  get  the 
following  estimates 

|e(3)|  .6(3)  ,   |Ae(3)|  <  gOJo  |iog  h|  . 

Ax     "" 
(2.1;5) 

|^e(3)|  <  5(3)0 Jiog  h|  . 

Ay  "     ° 

The  above  bounds  can  be  used  to  get  from  (2.6) 

(2.I16)         \A,y^K),^e^^^)\   <  a(^) 

where  5^^^  =  i;5+65AtN^  I  log  h|+I|.6(N^A^IloS  hI)^+5(AtN  Ilogh|)3» 
Prom  (2.i|6)  we  shall  get  the  bounds  for  e^'-^,  —e^^^,  — e^^^, 

A^         Ay 

Hence  vre  can  continue  this  process  a  finite  number  of  times. 
It  Is  obvious  that  after  n  time  steps  we  shall  get  the  follow- 
ing formula 

lAhe^^'^-k^e^^M  <  6[(;)  +  (^)lIoAt|log  hK(5)(N^At  • 

(2.1+7)   •  llog  hf)2+...  +  (^2i)(NQAt|log  hl)"-2+(^)(N^At  ' 
.  Iloghl)^-^] 


-i 


23 

or  f  r,  \ 

lr.\         O     (rs\  (d+N  At|lOg  hD'^-l 

(2.1^8)    lAh^^'^^-k^-^''^  I  <  S-{ 2f±l.^ . 

At  this  point  we  shall  use  the  Law  of  the  Mean,  i.e.,  consider 
N  At  I  log  h|  =  y,  hence  f(y)  =[(l+y)"l]/y  is  continuous  for 
y,<y  <  yp  and  differentiable  for  y,  <  y  <  ypj  therefore  we  can 
write  f(y)  =[(l+y)'i.l]/y  =  [(l+y)"-(l+0)^]/y  =  nd+Qy)"^"^  where 
0  <  e  <  1.  Hence  (2.L|.8y  becomes: 

lA.he^^^-k^e^^'M  <  6n[l+9N^At|log  hU'"-^  <  5n(l+N^At  I  log  h|)""^ 

.  n        N  Atnllog  h( 

"  l+N^AtUog  hi  l+N^AtUog  h| 

but  n/\t  =  T,  hence 

and."""""""')"    -"""  i°s  h 


lAh^'^'-k^eC^'l  < 


n     '      ■^.  5ne 


l+N^A^Ilog  h|    I-NqA*  log  h 

(2,1+9) 

-N  T 


"■^  '  "■  1-N^At  log  h 

Substituting  for  6  =  5^^  its  value  from  (2.8)  and  recalling 
that  A  =  A*A  and  T  =  nA*  we  shall  obtain  from  (2.i;9)  the 

following  inequality: 

-NT 

1-N^At  log  h 
or  ,^  1-N  T 

(M^+M^)Th   ° 


|At,e'">-k2o("'|  <1- 


N^Xh  log  h 


if  we  let  h  — >  0  and  1-N^T  >  0 

(2.50)  "T  <  r- 

o 


LJUlri  porUA.T- 


vUil.  j) 


rl  -  ^■■ 


''{Ii5  ROl;ct/\.T/r 


i-^A. 


V;    dn' 


icT   .1A 


■■J. •;';:« 
,  ?.^    rnnsF     .1  >   o  >  o 

1:6  >    |^^'''"o 


:u->    > 


I  (.'t), 


H ;..;.! 


-hr)n.3 


.--  >  P'^'o-'::-^'^' 


C/-J 


«*f- 


^"^-.Al 


^^|i.5) 


gni -1 -..::?  vx 


\ ) 


•   If. 


:^'?,s.) 


*Wfev 


2k 


we  obtain  that 

(2.51)        lAh^^''^-^^^^^''^  I  <  V^"'^°'^ 

hence  if  (2.50)  is  satisfied  the  Heliuholtzian  converges. 


3,  Estimate  for  the  Helmholtzian  (Initial  Error  E^°') 

In  the  previous  chapter  i-je  made  the  assiomption  that  the 
initial  error  of  the  Helmholtzian  is  zero,  i.e.,  /we ^ ^ ^ -k  e^°' 
=  0,  We  shall  novj  proceed  vjlth  a  convergence  proof  without  the 
above  assumption,  i.e.,  we  shall  consider  the  Initial  error  of 
the  Helmholtzian  to  be  given  by 

(3.1)  IAh-^'^-k^^^°^l  <E^°^  • 

From  (3.1)  by  previously  mentioned  and  in  part  II  in  de- 
tail developed  methods,  we  obtain 

(3.2)  |e^°M  <e(°^  ,   |^e(°M  <C  E^^Mloghl  , 

Ax 

l4-e(°)|  <C  E^^^log  h|  . 

Ay 

By  (2.6)  for  ra  =  0,  i.e.,  for  T  =  A^  ^®  have 

Ai,e(^^-k2e(l)  =  AtT(^).At[(f-e^°M^-(^(°h^] 

Ay        Ax   Ax        Ay 


(Ah*^°^-k2^^°))+[i+/;^t(^l^°h^ 

Z^       A 

(3.3) 

-  At(— l^°^)~3(Ah®^°^-k^®^°^) 
Ax        Ay 

A^e^^^-k^e^l)  =  AtT(^)-^K^°^L^°^ 


r    il^.(0)„    ^ 


r\v    re    •     •       -ion  Il^ris   6W     «0  « 

■l'i^J^.'.iil'\.1.0\i     yJTiw 


(o),. 


'.vl      f^.r\ 


•ol    fd.G) 


c^ . 


f'3)j.5.(0>^.v/.>;T'*   \    •■=    ^ 
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where 


Ay  ^ 

(3.5)  l(°)  =  u+At(^l(°))^-Z\t(^I^'^)^](Z:^.e(°lk2e(°)). 

Ay        A-        Ax        Ay 

Using  reasoning  similar  to  that  of  chapter  2  we  obtain  the 
following  upper  bounds  for  K^°'  and  L^°',  i.e., 

(3.6)  |K^°M    <  N^AtE^°Mlog  hi 

(3.7)  |L^°M  <max|Ahe(°^-k2e(°M  <  E^^^ 
with  the  following  restriction 

(3.6,  .  =  4i  < 1 


max 


I A^  Ay        J 


From  (3.3)  using  (3.6),  (3.7)  and  the  fact  that  A^I^^*)!  <  ^ 
we  have 

(3.9)  lAh^^^^-^^e^^^l    <  5+E(°HjQAt|log  hl+E^^^    =  E^^^ 

(3.10)  IAhe^'^-k'e(l)|    <e(1)    . 

From  (3.10)  by  previously  employed  methods  we  obtain: 

(3.11)  le^^M    <e(^^    ,      l^e^^M    <C^E(lMlogh|    , 

i4_e(l)|    <C   E^^Mloghl    . 

Ay 

In  a  similar  manner  as  before  from  (2.6)  for  m  =  1,  i.e., 
T  =  2 At  we  get 


»  '     * 


o ; 


A 


(o) 


Ui,c; 


{':.C} 


5  '-■-    , 


,P)   v^tnij   (C»n     m^'^ 

....  ■    ;  i 


:!  .1  i    .   I.,  i  -j 


i^i   -r>r! 


^03    ^Is 
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(3.12)  lAh^^^^-l^e^^M  <  At|T(^')|  +  lK^^M  +  |L^^^' 

where  this  time  K^  '    and  L   '  have  the  following  upper  bounds: 

(3.13)  IK^^M  <  N^AtE^^Mlog  h| 

and  consequently 

(3.15)   IAh®^^^-k2e(2)|  <  5+N^AtE^^Mlos  h|+E^l^  . 

Substituting   the   value   of  E^    '    in   terms   of   the   initial 
error  E^°^    from  (3,9)    into   (3.15)   we  obtain: 

lA^e^^^-k^e^^M    <  26+5N^AtUoS  h| +E(°U2E(°)N^At 
(3.16) 

.    Ilog  hl+E^CN^AMloS  h|)2   =  e(2) 

(3.17)  lAh^^^^-k^e^^M    <  e(2) 

and  from  the  above  by  methods  developed  in  part  II, 

(3.17)  |e(2)|    <  g(2)    .      |A.e(2),    <cE(2)|iogh|    ; 

Ax 

|^e(2)|    <  c   E^^^log  h|    . 

Ay  ° 

As   a  next   step  in   the   iteration  process  we  get 

lA^e^^^-k^e^^M    <  35+3SN^At|log  hl+eCN^AHlog  hl)^ 

(3.18)  +  E^°U3E^°^NQAt|log  h|+3E^°^(N^A* 

•  |log  h|)^-+E^°^(N  Attbg  h|)^  =  E^^^  . 

o 

It  is  not  difficult  to  see  that  continuing  the  above  itera- 
tion process  a  finite  number  of  times,  i.e.,  for  T  =  nA**  w® 
shall  get  the  following  formula  governing  the  error  of  the 


''(]- 


'^.^ 


;-o  > 


(►■'.fj  mo'sn   ^"'-'S  no-j/ra 


^^^o^r^^^a^ 


fa.c.c\ 


■r.uocfr.  £;■;;< 


/,r.^         A 


I      iv-.::.; 
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Helmholtzian: 

.|106  h|+(g)(N^At|log  h|)2+...  +  (^)(N^At  • 

•llog  h|)^] 
or 

(3.19)  lAh^^^-k^e^^h  <  6— ^^-^-^.^--^ W°^l+N  At  . 

N^AMlos  h| 

•  llog  h|)"  , 
The  first  part  of  (3.19)  which  \^e   r-hall  denote  by  P  is 
identical  with  the  right  side  of  (2.Ii.8)  and  its  estimate  is 
given  in  (2.Ij.9),  i.e., 

(1+N  At  llog  h  1)^-1     6nh'^°^ 

(3.20)  P  =  5 2=t <   onn 


NqA^IIoS  1^1     ~  1-NoAt  log  h 
We  shall  denote  the  second  part  of  (3.19)  by  S,  i.e., 

S  =  E(o>(l«„AtUog  hi,"  =  ,(o)„/oAtn|losh|^„ 

but  nAt  =  T,  and 

hence 

(3.21)  S  <  g(o)^  %'  , 

But 

iAhe^^'^-k^e^^M  <  P+S  . 

Therefore  by  (3.20)  and  (3.21)  we  obtain 


.'  •■'  v  f  J .. 


I'-- 


»  .-t  * 


5 


li;    /-i,;; 


H-£) 


.'  r:  ^     "  '\ 


n. 


rf/X^n  .j-.trcf 
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-^o'^     ,  .  -NT 


-NT 
0 


l-N^At  log  h 
Substituting  the  value  of  6  from  (2,8)  v/e  obtain 

-NT 
^r^^  o    i^\  (M  At+M,h)Th  °    ,. 

I-NqA*   log  h 
but  ^t  =  Ah,   hence  INT 

^  l-NQ;Vh  log  h  ^ 

^  l-N^-;Ui  log  h       ^ 

o 

If  E^°^  =  0(h)  and  T  <  1/N^  then  the  error  of  the  Helm- 

holtzian  — ->  0  as  h  — ->0,  and  the  growth  of  the  error  is  gover- 

ned  by  the  above  formula,  or 

(3.23)  IAhe'"'-k2e(T)|  <Ky-V 

Where  (T)  is  used  Instead  of  (n)  to  indicate  the  reached  time. 
By  methods  developed  in  part  II  we  shall  get 

(rp)  1-N-T 


W'^\    <Kh   ° 


ij. ,  Estimate  for  the  Helmholtzian  Extending  the  Total  Time  Be- 
yond Previous  Restrictions 
In  the  previous  chapter  our  formula  governing  the  error 

was  valid  under  the  restriction  that  total  time  T  <  l/N^  and 

—    o 

that  the  initial  error  E^°'  =  0(h).   Our  next  object  is  to  show 
that  it  is  possible  from  previous  results  to  develop  a  formula 


•'ir^- 


l^^^^o-^^^^" 


-^'JJjw 


\   3;*d 


(':^»i:) 


>  T 


-    ( 


•J.U     ^>-;.ii.  iu'J'.'l     JV'V'JfJ     €txi.-     7 


^■^Ti-^- 


••. ; 


O  c*0'O.L  v' '-''"-' i^'     K,  J 


;f>H  stroivev 


l&V    ^PA 


■i   ix>C/j 


j,ij'^  ■■ 
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which  would  show  us  that  our  proposed  finite  difference  scheme 
converges  even  for  a  larger  time  than  T  <  l/N  •  To  do  so  vje 
shall  proceed  as  follovjs.  Consider  that  we  reached  the  time 
T  <  1/N  which  from  now  on  will  be  denoted  by  T,  and  called  the 
first  "time  layer,"   In  the  previous  chapter  the  respective  er- 
ror estimates  for  the  Helmholtzian  and  the  error  itself  were 
given  by 

(T^)     1-NT,       p(o) 

where  h,  denotes  the  width  of  the  interval  ^x  =  /\y  =  h,  for 
the  first  "time  layer,"   The  above  formulas  are  valid  under  the 

restriction  that  T^  <  l/N^  and  E^°^  =  O(h^), 

(T  ) 
Let  us  consider  E^  1'  as  our  initial  error.  By  previously 

discussed  iterative  error  analysis  we  shall  obtain  the  follow- 
ing error  estimate  for  the  second  "time  layer"  with  interval 
size  /^x  =  /\y   =  h^ 

and 

(T,+To)      1-N  T^       T 
.^         Ul2|^v.   o2^prp4.  I 


(T,+To)      1-N  T^        J'^l'' 


(T,) 
Prom  the  above  formulas  it  is  obvious  that  E     has  to  be 

(T-,) 
of  order  hp,  i.e.,  E     =  ^^^2^  ^^'^   ^2  -  ■^^^^o'     ^^   *^^^  point 

we  shall  establish  a  relation  between  hp  and  h^  and  find  an  er- 
ror estimate  for  the  Helmholtzian  at  the  endpoint  of  the  second 
"time  layer"  in  terms  of  the  interval  h,.   It  is  not  difficult 


U.o   nSf.r 


•^'/Q 


.'  r.c^  ■"{'•' 


■3  bee 00':  ■% 


"T      ".-Tv    , 


■JB0    "IQ'I 


:..'.) 


(K,J) 


.1  ^  :■ 


■X  > 
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(T,) 
to  see  from  (li-.l)  that  for  E     to  be  of  order  hp,  taking  into 

account  that  E^°'  =  O(h^) ,  the  relation  i-ie   are  looking  for  is: 

1-N  T, 

ik.B)  h^  °  ^  =  h^  . 

Prom  (i}..3)  and  (l+.l)  we  obtain 

1-N  T, 
(T,+Tp)   p  (T^+Tp)         1-N  T,  h^   '^  "^      1-N„Tt 

E(o)h';^o'i 

or                                     l-N^T^  1-N  Tt 
(T,+Tp)   p  (T,+Tp)         1-N  Tp  CTT,hp   °  ^h,   °  "^ 
fAh^   ^  ^  -k^e   1  2  ,  <  ,^^^^^   o  2^  1  1  2  ^^ 1 . 


'2 

,  V  1-N  T^  l-N^T, 

By  (1;.5) 

1-N  Tp     (1-N  T,  )(1-1LT5)  • 
(i+.7)  hg  °  ^  =  h^^   °  =^     °  2'  ^ 

Substituting  (i+.5)  and  (i+.Y)  into  (14-. 6)  we  obtain 

and 

(11.9)  |e   1  2  ,  ^  g  2   ^ 

To  get  an  estimate  at  the  end  of  the  third  "time  layer"  we 

(Tp) 
consider  the  initial  error  to  be  E     and,  proceeding  in  a  si- 
milar manner  as  before,  we  obtain  the  following  estimate  for 
the  error  of  the  Helmholtzian: 

(T  +Tp+T  )     (T^+Tp+T.)      1-N  T-      „^^2^ 
(^.10)  I  Ah^  ^  ^  ^  -k^e  ^  2  3  I  <  ^3  °  \t^+  \—) 


■^'TO     ^O     Cf^     ■': 
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(•'^jjj)   trot'? 
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;^J)    c^n 
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where  h.  denotes  the  width  of  the  interval  /\x  =  /\y   =  h-  for 

(To)         ^ 
the  third  "tinie  layer,"   It  is  clear  that  E     has  to  be  of  or- 
der h-  and  T.  <  l/N  .   or  this  purpose  using  (i|..8)  and  (k*$) 
we  establish  the  following  relations  betvjeen  h^  and  hp,  and  be- 
tween h^  and  h, ,  respectively: 

I-NT5    (1-N  T-,)(1-?L T^) 
(1^.11)         ^3  =  i^      =  ^1 

and 

1-W  T,     (1-M  T, ) (1-N  T3) (1-N  T J 
(1^.12)       h^  °  3  =  h^   01     0  2     03. 

Using  (i;.8)  and  ([;.ll),  (I|,10)  becomes: 

(l^•l3)  „(o) 

and 

(T,+Tp+TO      (1-N  T^  )  (1-N  T  )  (1+N  T.) 
(i^.ll}.)  |e^^2|<h3^   °-^     °2     03  [c^(T3_+T2+T^)  + 


h. 


-]  . 


If  we  continue  the  above  process  a  finite  number  of  times, 
suppose  k  times,  we  shall  obtain  the  follovring  error  estimate 
for  the  Helmholtzian  at  the  end  of  the  k*^  "time  layer" 
(T3^+T2+...+T  )     (V^2'■•••■'\^ 

(1-N  T-,)(l-N  Tp)...(l-N  T,) 
(U.15)  <  hi     -^     °  ^        °  ^ 

with  the  restriction  that  for  each  "time  layer"  T.  <  l/N  and 
E^°^  =  0(hi).  From  (I|..l5) 
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,4.16,    i^VV-*-.),  .   '^-Vi)<^-V2)-<i-W 

r.(0) 

.    [(T^+...+T^)C^+£^]    . 

If  we   take   the   restrictions   on  T   for  each  "time  layer" 

T      <  -=- 

1  -  \ 

T      <  — 

2  -  No 

• 

k  -  r 

o 

and   add   the   Inequalities,   we   get 

W...fTj^   <^   • 

From  the  above  inequality  we  see  that  it  v;ould  be  conven- 
ient to  choose  equal  time  steps  for  each  layer,  i.e.,  T^  =  T2  - 
=  ...  =  Ti,  =  T,  hence,  if  we  denote  the  total  time  for  all 
"time  layers"  by  T,  we  obtain  T  =  kT  and  ()4.l5)  becomes 


C^)      ?  fT'l,     (1-N  T)^  ^   g(o) 


but 


hence 


and 


where 


^(i-V)\,a-(™,A))^,^e-™< 


IAj,em-lc2e(T)|  <hf    (tV^) 


^vj 


o     1 


If  we  choose  E^°^  =  O(h^)  then  as  h^  — >0,  |e^^'|  — >0  (be. 


no 
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cause  e     is  a  positive  number  and  h-,     — ■>0  as  h,  — ^>  0)  • 

^   k 

Hence  our  scheme  is  convergent  for  T  <  ||-. 

0 


5.  Round-Off  Error 

In  the  previous  error  analysis  we  assumed  that  we  can 
solve  the  finite  difference  equation  (2.1)  with  infinite  preci- 
sion, i.e.,  the  effect  of  the  round-off  has  been  ignored.   In 
practice,  however,  we  compute  the  solution  of  (2.1)  rounded-off 
to  a  certain  number  of  decimal  places.   If  we  keep  the  number 
of  decimal  places  fixed  and  decrease  the  mesh  v;ldth  vje  could 
not  expect  convergence.  VJe  shall  show  that  the  scheme  is  con- 
vergent if  the  round-off  is  of  order  O(h^).   Let  us  denote  by  U 
the  quantity  U  after  round-off,  i.e.,  U  =  U+0(h^) . 

It  is  not  difficult  to  see  that  in  solving  (2.1)  for  the 
Helmholtzian  we  must  require  that  the  round-off  be  of  order 
0(h2),  i.e.,  j:/'^'^^   =^"''^^+0(h^)  where  /-/  =  Aur^^i*  Conse- 
quently we  establish  that  the  round-off  for  j|j  itself  should  be 
of  order  O(h^),  i.e.,  f^^   =  f^^+0(h^*-). 
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Part  II 

6.  Estimate  for  the  Error  e  =  i]/  -  ^ 

To  find  an  estimate  for  e^^  from  (2.7)  we  shall  employ  a 
special  method  given  in  [?]  v/hich  involves  the  construction  of 
an  auxiliary  function  and  is  called  the  method  of  majorants. 
However,  before  proceeding  with  the  above  mentioned  method  we 
first  need  to  prove  a  few  theorems. 

Theorem  1.   If  quantities  v,  ^at  all  interior  points  of  the  grid 
domain  G,  satisfy:  /\,  v-k  v  <  0  and  at  the  boundary  points 
V  >  0,  then  v  >  0  at  all  interior  points, 

p 
Proof.  Considering  /^^,  v~k  v  defined  as  before  by  (1.9)  we  as- 
sume the  contrary,  i.e.,  that  at  some  points  in  the  interior  of 
Gj^  the  function  v  assumes  negative  values.   Since  v  is  non-nega- 
tive on  the  boundary  there  would  be  an  interior  point  at  which 
v  would  assume  its  negative  minimum,  i.e.,  v  =  a  <  0  and  at 
least  one  neighboring  point  at  which  v  >  a,  i.e., 

^ki^^k-M,i'. 

^k^^^k-i,v: 
^k^^^k,v+i 

^k^:2^k,^l  ; 
then  it  is  obvious  that  in  at  least  one  of  the  inequalities  the 
equality  sign  does  not  hold,  hence  after  adding  these  inequali- 
ties and  taking  the  above  into  consideration  we  obtain 

^k+i,.e-'^k-i,|-'^k,j?+i^^k,i'-i-K£  >  ° 

and 

-h^k^v^^>   0 


iO 


ti-ig   ciiJ    iw   ■■^:-- 


^  r,   O         C^» 


V- 


• '   ;  \  J    -'ir  * 


a^ 


iVi;J 


Q    ax 


i  .-r,. 


:;'"|i'^^f""'i 


35 
because  by  assumption  v,g  <   0  and  adding  once  more  we  get 

which  contradicts  our  assumption  that  A,  v-k  v  <  0.  Hence  we 
see  that  our  assiimption  about  v  taking  negative  values  in  the 
interior,  which  lead  us  to  this  contradiction,  must  be  false, 
and  therefore  we  conclude  that  v  >  0  at  all  interior  points  of 

Theorem  2«  If  quantities  v,  ^  and  V,  ..  at  all  interior  points  of 
the  grid  domain  G,  satisfy  the  follox^ring  inequality:  A,^V-k  V  < 
<  ~|/\^v-k  v|  and  at  the  boundary  points  V^.  >  lv^.',|  then  Y^^  > 
>   |vt  ^1  at  all  interior  points  of  the  grid  domain  Gy^. 

Proof.  The  above  theorem  follows  directly  from  Theorem  1  if  we 
observe  that  the  inequality 

^j^v-k^v  <  -IAh^-^-^^v| 

is  equivalent  to 

A^^V-k^V  <  Ah^-k^v  and  Ah^-k^V  <  -(Aj^v-k^v) 
or 

(6.1)  Ah^^-^^-^^^^-^^  ^  °    ^"^   Ah^^'^^)-^^^^'*"^^  ^  ° 

and  that  on  the  boundary 

(6.2)  v^y-^kk.>o   ^^f.-^^ki^o- 

Considering  (6.1)  and  (6,2)  v;e  apply  Theorem  1  to  quantities 
V+v  and  V-v  and  conclude  that  at  all  Interior  points  of  G, 

V+v  >  0  and  V-v  >  0 
or 

(6.3)  |v|  <  V  . 
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To  estimate  the  error  e'^  from  |/\,  e^^-k  e^  '|  <  5  we 
shall  employ  the  above  mentioned  method  of  majorants  vjhose 
main  idea  is  to  construct  a  bounded  function  z(x,y)  which  in  G, 
would  satisfy  the  following  inequality  -(/\,  z-k  z)    >  6   and  get 
the  estimate  applying  theorem  2. 

We  consider  our  grid  doraain  G,  to  be  a  square  |x|  <  1, 

|y|  <  1  and  circumscribe  around  it  a  circle  of  radius  r  =  /2 

and  center  (0,0)  .  V/e  construct  the  following  fimction 

2   2 
(6.k)  z(x,y)  =  p[l-|--^] 

and  apply  the  operator  ( /\} -k  ),  i.e., 

(6.l}.a)  Ah^-I^^z  =  H(z)-R^^y^(z) 

where  H(z)  and  ^y^A^)    are  defined  by  (1.11)  and  (1.12). 

H(z)  =  z^^+z   -k^z 
XX   yy 

2   2  2   2 

H(z)  =  -(3-p-k2p[l-|--:|.]  =  -2P-Pk2[l-^-^]  . 

But  R,  ^  in  our  case  reduces  to  zero,  because  it  involves  deri- 
vatives of  order  higher  than  2,  i.e.,  ^1,0(2)  =  0,  Hence  (6.l4.a) 
becomes; 

AhZ-k^2  =  -2(3[l+^(l-^-^)] 

2    2   2 
-(/y^z-k^z)  =  2[3[l+|-(i.?^-^2-)] 

P 

but  the  factor  of  k  /2  never  becomes  bigger  than  one  and  it  Is 

>  0  in  G,  ,  hence  we  conclude 

-(/^.^z-k^z)  >  2(3 

and  to  achieve  that  -(^^  z-k  z)  >  5  it  is  enough  to  choose  p  to 
be  equal  to  5/2,  i.e.,  p  =  5/2,  and 
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(6.5)  -(Al,z-1^^2)  >  6  ; 

combining  (2.7)  and  (6,5)  we  obtain 

-(Ah^-l^^z)  >  5  >  lAh^^^^-k^e^^M 

at  the  interior  points  of  G,  .   On  the  boundary  z(x,y)  =  0  and 
if  we  assume  that  e  >  0  on  the  boundary  vje  are  able  to  apply 
Theorem  2  and  conclude  that  at  all  interior  points 

le^^M  <z 
but 


2   2 
2'   2   2^-2 
and 


z=£[i.^:f.-C.]  <5.  <  6 


(6.6)  |e^^^ I  <  5  . 


7 •  Simplification  of  the  Auxiliary  Problem  and  "Discrete  Po- 
tential Equation" 
In  this  chapter  vje  shall  be  concerned  with  the  solution  of 
the  following  problem 

,2^  _ 


(7.1)  ^  ^  0 


on  the  boundary  of  G, 
where  |f (P) |  <  6. 

Let  us  reduce  the  above  problem  as  follows: 
/\^e  =■■   k2e+f(P)  =  g(P) 

by  (6.6)  |e|  <  5  hence  |g(P) |  <  C-S  where  c  =  k  +1,  and  we 
shall  have  the  following  problem 

(7.2)  ^  "■ 

e  =  0       on  the  boundary  . 
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We  shall  show  that 

e  =  e^+e^ 

Is  a  solution  of  the  above  problem  if  e,  is  a  particular  solu- 
tion of 

and  ep  is  a  solution  of  the  follo^^^ing  boundary  value  problem 

Ah^2  =  0   V  eG^ 

(7  5) 

'  Op  =  -e,  on  the  boundary  • 

If  u(P,Q)  is  a  function  such  that 

(7.6)  Ah^(P.^)  = 

then  a  particular  solution  of  (7»^)  can  be  written  as 

(7.7)  e^^CP)  =  h^  2I2Z   u(P,Q)g((i) 

where  the  sum  is  taken  over  all  grid  points  Q,  vxlthin  the  boun- 
dary of  G,  .  For  if  we  operate  Ah  ®^  defined  by  (1.3)  on  both 
sides  of  (7.7)  with  respect  to  the  coordinates  of  P,  and  consi- 
der (7»6),  we  obtain 

(7.8)  Ah«i(P)  =  h^  ^  Z  AhU(P»^)g(<^)  =  s(P)    . 

Q. 

For  further  analysis  we  shall  need  the  estimates  of  (7.7) 
and  its  difference  quotients,  i.e., 

(7.9)  |e3_(P)  I  <  max|g(Q)  \h^   ^  ^I  |u(P,Q)  I 

(7.10)  l~e,  (P)  I  <  max|g(Q)  \h^   Z  Z.   l^u(P,Q)  | 

Ax  ^  A^ 

(7.11)  |^e,(P)|    <raaxfg(Q)|h2  ^l  Y.   |^u(P,Q)  |    . 
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To  obtain  the  above  estimates  v/e  actually  have  to  estimate  the 
following  sums 

(7.12)  h^  21  Z  |u(P,Q)I 

(7.13)  h^  ^21  'X;''^^'^^' 
il.lk)  ^^  HH   l~-u(P,Q)|  . 

Ay 

To  do  this  we  first  have  to  obtain  tho  bounds  for  the  so- 
lution of  (7.6)  and  its  difference  quotients, 

A,  Stohr  in  reference  [6],  B,  von  der  Pohl  in  [53  and  S. 
L.  Sobolev  in  [8]  gave  an  extensive  study  of  the  solution  of 
the  following  equation 

(1  X  =  r  =  0 
itherwise 

where  x,y  are  positive  or  negative  integers  Including  zero,  and 
the  symmetrical  difference  opero.tor  \/   is  defined  by 

(7.16)  \/vi{'x,y)    =   u(x+i;y)+u(x-l/y)+u(x,y+l)+u(x,y-l)-l;u(x,y)  . 

We  shall  use  soine  of  the  results  of  B,  von  der  Pohl  [53 

u 

and  A,  Stohr  [6]  to  obtain  the  estimates  for  the  above  mention- 
ed sums  by  showing  an  analogy  between  the  solutions  of  (7.6) 
and  (7.15). 

B.  von  der  Pohl  in  reference  [$]   verifies  that  the  follow- 
ing function  is  a  solution  of  (7.15) 

2ti    Sir       o  .^/ 

(7.17)  u(xVy')  =  fe  ^  4^  (  d^.-^MlTZl).  . 

Vie  shall  now  verify   that  , 

(7.18,        u,.„;,,r,,  =  ^  j  a+lj]  a*!—-^-^ 


(7.15)  Vu(x,y)  =  i 

I    0      oi 


/>_. 
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is  a  solution  of  (7.6)  ,   Let  us  apply  At-  on  (7«l8) 

(7.19)    Ah^  =  ^2  ^* ;  ^4    — 


O    O 


sin  <j)+  sin  vj; 


It  can  be  easily  verified  that 


f 


A 


(7.20) 


Ah^sin2[(x-?)J+(y-7/)|](  =  -^(sin2(},+  sin^^l,) 


h' 


COS  2[(x-5)~^+(y-t|)f]  . 


If  vje   substitute    (7.20)    into    (7.19)   w©  obtain: 

271     2ii  A  • 

(7.21)  Aj^u  =  -i^j  d4  j  dt-~cos  2[(x-^)i+(y-;.)|]    , 

0  0 

If  P  s  (^,   i.e.,    (x,y)    =   (^jYj),    (7.21)   becomes 

2Tt       271 


''  -  ;;?>'  =  ? 


hence 


AhU(P,^)    =  -^     for     P  =  Q   , 


and   if  P  ?^  Q  fro.m   (7.21)   we  get 

2ii     271 


A 


4 


(7.22)Ah^  =  "T  "T  !  "^M  dilr[cos(x-?)Jcos(y-rj)|-sin(x-5)jjsln(y-r^)|] 

^^  O  0 

but  i-jQ   consider  only  the  values  at  grid  points,  hence 

X  =  xh 


(7.23) 

and  (7.22)  becomes 


^   and     ^ 
y  =  yh       T/  =  ^/h 
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Ah^  =  -^'o)  \  cos[{X'%)^]A^  \  cos[(y-'\;)\!/]d\|/-  (  sln[  (x-?)(l))d(i). 


2% 

•    1  sin[  (y-'T))i{']d\i/ 


1  =  J^  i  sin(x-g;)<l' 
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hence 


"■  l;Tch^  x-^     x-5 

hence  Av.u  =  0  for   P  7^  Q,.     V/e  have   shown   that    (7,l8)    is   a  solu- 
tuion  of    (7.6)    and   if   into   its  i-'ight   side  we   substitute   for    (x,y) 
and    (C»'^()    its   values   from   (7.23)    we   obtain   for   the   solution  of 
(7.6)    the  following  expression: 

-   '  2ti     2ti  p     ,..  „,,         _  _^ 

'  "^^l       i      '^  sin'^d)+  sin'^ij/ 

0  0  T  T 

Before  proceeding  with  tha  estimates  v;e  shall  establish  a 
few  properties  of  the  above  solution. 
Prom  (7.18)  it  is  obvious  that 

(7.25)  u(x,y;^,-y)  =  u(^,?);x,y)  . 
Comparing  (7.17)  and  (7.2li.)  vje  see  that 

(7.26)  u(x,y;0,0)  =  u(xVy;0,0)  =  u(x',y)  . 
From  (7.2lj.)  it  is  obvious  that 

(7.27)  u(0, 0,0,0)  =  0  . 
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Q,     A  Bound  for  the  Non-HomoReneoiis  Case 

Next  we  seek  an  estimate  for  (7.12).   It  is  obvious  that 
considering  (7.25)  and  (7.26)  \<ie   can  get  it  by  estimating  the 
following  equivalent  sura: 
(8.1)  ^^  HI.    u(x,y)  . 

The  way  is  now  prepared  to  employ  the  following  bound  for  u(x,y) 
given  by  A.  Stohr  in  reference  [53: 


M 


(8.2)  |u(x;yV^log  s/5^-^:bs  2-2^01  <  ^2-^  , 

X  Ty 

Where  x  +y  7^  0,  C  -  Euler's  constant,  and  M  -  a  positive  nura- 

ber  independent  of  x  and  y. 

From   (8.2)   vje   obtain 


(8.3)  |u(x,y)|    <   [log /x^+y^|+K+,-r2-^%^ 

X   +y 

where  K  =  j^-log  2+™^C. 

Having  in  mind  (7,27),  we  should  sum  (8,3)  over  the  follow- 

ing  region  1  <  x  +y  =  r  <  R  .  Hovrever,  let  us  subdivide  the 

above  region  in  two  parts  as  follovjs: 

(8.L|.)  1  <  r^  <  (R')^ 

(8.5)  (R')^  <  r^  <  (R)^ 

and   take   the   following  bounds   for  u(x,y),   respectively: 

(8.6)  |u(x,y)|    <  ^|log  yP+y^l+K+M     for      1   <  r^   <   (R')^ 

(8^7)       |u(x",y)  I    <  ^|log  v^?T^|+K+-ii.^  for    (R*  )^   <^v^   <  (rV    . 

'^'^  X   +y 

Using  (8.6)  and  (8.7)  vje  obtain  the  following  estimate  for 

(8.1) 
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h2  Z  Z  Iu(^~  S)  I  5  h^  Z  Z  tilllos  ?I*L1  + 
(8.6)  ^^^«* 

where  L  =  K+M, 

Returning  to  the  original  plane,  i.e.,  taking  r  =  r/h,  we 
obtain  the  following  subdivisions  corresponding  to  (S.li.)  and 
(8,5),  respectively: 

h  <  r  <  R« 

R»  <  r  <  R 
and  after  choosing  R»  =  1  our  range  becomes 

h  <  r  <  1 

1  <  r  <  R   . 
Hence  we  obtain  the  following  expression  for   (8,8). 

h2  Z  Z  Iu(x,y)  I   <  h2  21  r  [^llog  |1+L]   + 

h<r<l 

*h2ZZ[^|loggl«-^i 
l<r<R    '^^  v^ 

h^  2:Z  lu(^,y)I  <h'^T.T,  liog  ^I-  lih^  log  hUZ  1  + 

^^         h<r<l  ^^  l^<r<l 

+  Lh2  X  Z  1+  ^  h^  ^  £:   [log  r|-  ^  h^  log  h. 

h<r<l         '^^         1^<^ 

2 

•  Z  Z  l+Kl-i2  ^  X  l+h^M  Z  Z  ^ 
l<r<R  l<r<R  l<r<R  r 


h^  ZZ   |u(x,y)l    <  ^^h'^  HZ   |los  ^1+  ^  h^UoS  ^(^111+ 

'^^         h<r<i  '^^  h<r<R 


2 

(8.10)  +  Lh^  Z  Z  1+i^^  Z  Z  i+ti^M  z  z  \  + 

h:^<i  1<^<R  ^S"^^  ^ 

+  ^h2  HZ   Hog  ^1    • 
'^^  l<r<R 
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We  shall  first  approximate  the  following  sura 

(8.11)  h^  21'  ZI  I  log  r  I  . 

h<2?<l 

Let  us  now  consider  an  h  x  h  square  in  the  subregion  h  <  r 
<  1  of  our  grid  domain  G,  ,  and  let  us  associate  with  each  square 
a  mlnlma'Ti  value  of  |log  r|  at  a  grid  point.   It  is  obvious  that 
we  have  to  take  the  farthest  point  from  the  origin,  i,e,,  in 
the  first  quadrant  the  point  at  the  upper  right  hand  corner  of 
each  h  X  h  square;  in  the  second  quadrant  --  the  point  at  the 
upper  left  hand  corner;  in  the  third  quadrant  —  the  point  at 
the  lower  left  hand  corner,  and  in  the  fourth  quandrant  the 
point  at  the  lower  right  hand  corner,  as  indicated  in  Fig,  1, 

Then  it  is  obvious  that 

jl  I  log  r|dxdy  >min|log  rjh   • 

h 

Siairaning  all  the  h  x  h  squares  of  the  region  h  <  r  <  1  we 
obtain 

(8.12)  Z  Z  f  Uos  r|dxdy  >  h^  Z  Z  ilog  r|  . 

The  sum  on  the  right  hand  side  is  the  sum  over  all  grid 
points  except  the  points  on  the  horizontal  and  vertical  axes. 
Hence,  to  make  the  right-hand  side  of  (8,12)  the  summation  over 
all  net  points,  we  add  to  both  sides  of  the  above  inequalities 
the  following  expression  ^h  log  h+b  h   (where  I|./h  =  number  of 
points  on  the  vertical  and  horizontal  lines  and  b  a  positive 
constant) . 

J  J  I  log  r  I  dxdy+^h^log  h+b^h^  >  ^^  Z  Z  I  log  r  | 
h<r<l  ^^^  ^'P* 
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lim  [  I  (39  ( I  log  r|rclr+l4h  log  h+b  h^   >  h^  ^  Z]      |log  r| 
ti'>0    I       I  °        ~  all  n.p. 

I  d©    I  log  r|rdr  >  h*^  21  !^     |l°S  r*! 
1       I  all  n.p. 

(6.13) 


o   o 


2n      1 

h^  2]  21   tlog  r|  <  i  d© 1 1  log  r|rdr   , 
all  n.p.         >       ' 


o       o 

Now  we   shall  estlraate   the  remaining   suras    in   (8.10) 


^h^ilog  h|   ZZ  1  -  ^,,-^^llog  h|[if.|+f  ] 


(8.110 


l<r<R 


=  ^|log  h|+fh|log  h| 


(8.15)  pW[log  h|   Z  Z  1  =  ^Ipllog  h|+d|h|log  h| 

(8.16)  Kh^  T;  Z  1  =  Kh^-[-'^+-^]   =  bn+^;h 


27t      R 


(8.17) 


l<r<R  r'^  o       1 


271      R 


(8.18)  h^  Z  Z  |los  iH  '^     5  ^^('^"S  r|r  dr 


l<r<R 


o       1 


Hence  by   (8.13),    (8.15),    (8.16),    (8.1?)    and    (8.l8)   we  obtain 


for    (8.10) 


V 


271     R 
f 


|u(x,y)|    <       dojilog  r|rdr)  +  d2|logh|+d2h|log  h|+b'£h+ 


) 
o        0 


+b>d^h^ 
0      o 


li^  Z  Z   Iu(x.y)l    <  27iR^[i2|_?.-l]+b'^Vd2|log  hl+d^hllog  h|  + 


+b!^h+d  h^. 
1       o 


(8.19) 

h^  ZZI'^(-"^»y)l    <  <3^+d2|log  h|+d|h|log  h|+b£h+d^h2    . 


:./]    ^ 


^rrl'^b  '  ]  nil 
1-    ()'<--' 


c-  v.a 


!■    1 


-■r.  ',) 


\A  vi. 


:?<"'-h- 


icx.ii) 


' .  ^ 


'-jj-  <; 


d 


Now  we  are  read'y  to  obtain  a  bound  for  e,(P),  By  (7.9) 

(8.20)       |e^(P)|  <  max|g((^)|h2  Z  Z  l^(P»'^)l  • 
But 

S(P)  =  k2e(P)+f(P) 
hence 

|g(P)|    <  k2|e(P)|+|f(P)|    . 
By    (1.6)    and   the   fact   that    if(P)|    <  5  we   obtain 

is(P)  I    <  °o^  vjhere   c^  =  k^+1    . 

Therefore  considering  the  above  bound  for  g(P)  and  (8,19),  from 

(8.20)  we  get 

|e3_(P)l  <  CQ5[d^+d2|log  hl+djhllog  h| +b'^^h+d^h2]  . 
Since  by  (2.8) 

6  =  MQ(/;xt)^+M^^th+0(h^) 
we  shall  have 

(8.21)  |e^(P)|  <  c^Sllog  hl+0(h2)  . 


9«  Another  Bound  for  the  Non-Homogeneous  Case 

In  this  chapter  vje  shall  find  an  estimate  for  (7,13).  By 
a  similar  argument  as  in  chapter  8,  it  is  obvious  that  (7.13) 
can  be  approximated  by  estimating  the  following  equivalent  sum 

(9.1)  ^^  E  Z  E  ^«(^'»y) 

since  it  can  be  easily  seen  from  (7.17)  and  (7.21].)  that 

(9.2)  ^-  u(x,y)  =  i  ^^^u(x,y)  . 

/Xx  1^  A^x 

II 

A.  Stohr's  paper  [6]  gives  an  estimate  for  the  above-mentioned 
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function,  i.e., 

|u(x+l,y)-u(x,y)  |  <  ^jlog  ^(x+D^+y^-log  v/x^+y^|  + 


(9.3) 


^p  -x,p    /vp 

for  x'^+y'^  =  r*^  >  1. 


M        M 
,'-■  »2  "2   "2  -2 


We  shall  subdivide  our  region  in  the  following  way: 

2  <  r  <  (H' ) 

"^    p      rv  p        ^   P 

(R»)'=^  <v'^<   (R)^  . 


'V.-P   '^'P 

Since  x  +y  >  1,  it  is  obvious  that 


M   ^. M ^    M     M 


and  we  shall  take  the  following  estimates  for  different  regions; 

|u(x+l,y)-u(x,y)  |  <  ^|log  \/(x''i-l)2+y2-iog  v/x^+y^|  + 

^■^•^^  ^  __M M 

Vx^+y^   y(x+l)2+y2 
for  2  <  r^  <  (R')^  and 

|u(x+l,y)-u{x",y)|  <  ^|log  /(x+D^+'J^-log  yx^+y2|  + 
(9.5) 


_j_  M     ,   M 

(x+] 
p   ^p    '>^  p 


(x+l)2+52  -2^-2 


A 
Considering  (9.2)  as  I'/ell  as  (9.i+)  we  obtain  for  ^=*-u(x,y) 


for  (R')   <  r*""  <  (R)^, 
Considering  (9.2) 
the  following  estimate  in  2  <  r  <  (R') 

|^u(x,y)|  =  i|u(x;i,y|-u(x,y)|  <  1  1  |log(l+§f:^i)  |  + 

M  _^      M 


7x2 +'y^   y(x+l)2+y2 
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at 


!  f  •• 


Introducing  polar  coordinates,  i.e.,  letting 


X  =  r  cos  a 
y  =  r  sin  a 


we  get 

t^  „/..  „\  I   ^  J.  i   J-   t^f  OPS 

Ax  I       "^       "^      V^-+2r  cos  a+1 


1|  1  ar'cos  a+l^M_^ M_ 


'Z 


but 


hence 


<  i  i  1  2IH-1  ,  M  ,     M     '] 
I.    ^    ^   yr2-2r+l  ) 

I^ufx^n  <  1^2^.1+ 1.2.+  M  +  JL.1 

|—     uix,y;  I    <  T-  \2'n~li^^~'^-       f 
^x  ^  (  ^^t  r     i*-^  r'^     r     r-l  J 

1.1 


Il-Ttr     ""  linr' 


for         r  >  1     , 


I A   ..,_.   __v  ,     ^  Ir^o  .     H 


/\x  r      r-l 


• 


Returning  to  the  original  plane,  i.e.,  letting  r  =  rh  we  obtain 

(9.6)  lf-u(x,y)i  <-l^^ 

in  2h  <  r  <  1  if  vje  choose  R'  =  1,  As  a  next  step  we  shall  es- 
timate 

(9.7)  h2  Z;  Z  l^u(x,y)  I  <  h2  X:  Z  ^+^2  T  T  ^     ' - 

2h<r<l  /\j:  2h<i^l^     2h<r<r  ^ 

We  shall  first  approximate  the  following  sum: 
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As  before  we  consider  an  h  x  h  square  in  the  region  2h  <  r 


,  2  r-  Y-  "^o 


<  1  and  let  associate  with  each  square  a  minimum  value  of  Pq/i** 
It  is  obvious  that  we  have  to  take  the  farthest  point  from  the 
origin  and  exactly  as  in  chapter  8  we  will  obtain  for  every 
h  X  h  square  the  following  inequality: 


jj  ^dxdy  >rain(^).h2  . 

9 


h 


Adding  the  above  inequalities  for  all  h  x  h  squares  in  the 
region  2h  <  r  <  1  we  obtain 


ij 


Po...._  _  .2  ^  ^  Po 


ZZ  jJ^dxd,>h^ZZ 

D 

where  the  sum  on  the  right  hand  side  is  over  all  the  grid 
points  of  the  region  2h  <  r  <  1  except  the  points  on  the  hori- 
zontal and  vertical  lines.  Hence  if  we  add  to  both  sides  the 
(max  value  of  P_/r)  x  h  x  (number  of  net  points  on  vertical 
and  horizontal  lines)  we  obtain: 

J   J  ^o     2h^^h    -   all  n  pts.^ 


2h 


.2 P' 


2Tip„(l-2h)+2(p^-p  h)  >  h-^^  Z  Z  -jr 

all  n  pts. 


or 


2  ^  ^-  Po 


(9.9)  h^  Z  Z  ^  <  '^o^V  • 

all  n.p. 

As  a  next  step  we  shall  evaluate 
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(9.10) 


h^ z z  A' 


As  before  we  associate  v/ith  each  h  x  h  square  a  minimum 
value  of  M/(r-h)  vjhlch  v/ill  occur  as  before  at  the  farthest 
point  from  the  origin,  and  as  previously  we  v;ill  obtain 


If 
h 


^dxdy  >  mini^)   x  h^  ; 


again  summing  all  such  inequalities  for  all  the  squares  we  ob- 
tain 


2:^1 


A'^^y  >  h2  z  z  -^ 


but 


g^dxdy 


2h<r<l 


^.ZZlf^-a, 


hence 
(9.11) 


f  (  !2z!l .  h2  ^  T-  M 

J  J   r/2   ^  ^  ^  2-  ?IK  • 
2h<r<l 

The  sum  on  the  right  of  (9.11)  is  over  all  the  grid  points 
except  the  ones  on  the  horizontal  and  vertical  lines. 

Adding  to  both  sides  of  (9.11)  the  following  expression: 
h  X  (no.  of  points  on  horizontal  and  vertical  lines)  x  max  x 


M 


(pih^»  we  have 


2v.      1 

\ 

o       2h 


2|a«jg„ar<,i.i,l,>h2rZ   5^ 


-    ,    M 


h''  ZZ  iA  <<+<^ 


(9.12) 

Combining    (9.9),    (9.12)    and    (9.7)   we   obtain 
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(9.13)         ^^  lET.   I— u(x,y)|  <m,+n^h 

2h<i<l  ^x  ^     ^ 


where  n,,n,  are  constants. 

However,  we  need  an  estimate  of  (9.2)  in  the  whole  domain 
0  <  r  <  1,  therefore  we  shall  next  find  an  estimate  for  the  re- 
maining part:  0  <  r  <  2h,  VJe  estimate  the  difference  quotients 

II 

from  numerical  values  of  the  function  given  by  A,  Stohr  as  well 

as  by  B,  van  der  Pohl  in  [6]  and  [5],  respectively.  Hence  we 
consider  the  following  sum 

h^  ZZ  l~-u(x,y)  I    =  h2  3:  X   Au(x,y)  |  + 
0<r<l     Ax  0<r<2h  Ax 

(9.11+)  ~-      -^  --      _. 

+  ^2  2:E   l4^u(x,y)|    . 
2h<r<l  ^x 

The  estimate  for  the  second  siim  on  the  right  hand  side  is 

given  by  (9.13).   The  first  sum  will  consist  of  the  following 

quotients  in  I^*,  II"*^,  III^^  and  IV*^  quadrants,  respectively 

(see  Pig,  1) • 

(9.1^)        |[u(2h,0)-u(h,0)]  =  i[u(2,0)-u(l,0)] 

(9.1^,)       |[u(2h,h)-u(h,h)j  =  i[u(2,l)-u(l,l)] 

(9.1^.)        i[u(h,0)-u(0,0)j  =  i[u(l,0)-u(0,0)] 

(9.1^j)        |[u(h,h)-u(0,h)]  =  i[u(l,l)-u(0,l)j 

(9.IIg^)  5[u(-2h,0)-u(-h,0)]  =  i[u(-2,0)-u(-l,0)] 

(9.11^)  i[u(-2h,h)-u(-h,h)]  =  |[u(-2,l)-u(-l,l)] 

(9.IIq)  i[u(-h,0)-u(0,0)]  =  i[u(-l,0)-u(0,0)] 

(9.11^)  ^[u(-h,h)-u(0,h)j  =  i[u(-l,l)-u(0,l)] 
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(9.111^)  i[u(-2h,-h)-u(-h,-h)]   =  ^[u(-2,-l)-u(-l,-l)] 

(9.111^)  |[u(-h,-h}-u(0,-h)]  =  i[u(-l,-l)-u(0,-l)] 

{9.IV^)  |[u(2h,-h)-u(h,-h)]   =  |[u(2,-l)-u(l,-l)] 

(9.IV^)  ^[u(h,-h)-u(0,-h) j   =  |[u(l,-l)-u(0,-l)]    , 

A,  Stohr'in  [6]   gives   the  following  relations   for  u(x,y): 

u(x,y)    =  u(x,-y)    =  uC-x/y) -u(-x,-y)    =  u(y,x)    =  u(-y,x) 

(9.15)  ..    ,..  ^    ^'  ■ 

=  u(-y,-x)   =  u(y,-x)  • 

Due  to  which  the  difference  quotients  only  in  the  first 
quadrant  have  to  be  estimated,  since 
(9.1^)  =  (9.IIg^) 

(9.1^)  =  (9.11^)  =  (9.111^)  =  (9.IV^) 
(9.1c)  =  (9.IIc) 

(9.1^)  =  (9.11^)  =  (9.1X1^)  =  (9.IV^)  . 

It 

From  A,  Stohr's  table  in  [6]  v;e  have 

(9.16)  i[u(2,0)-u(l,0)3  =  ^[1-|-J]   =  ^ 

(9.17)         i[u(2,i)-u(i,i)]  =i[.i+|-ij  =:^ 

(9.18)  i[u(l,0).u(0,0)]  =  i[|-0]   =  ^ 

(9.19)  i  [u(l,l)-u(0,l)]  =  E^^'l^  =  ^  . 

Hence 

?  A  5i2b-,+2b,+l+b-+Ub.    1 

(9.20)  h^Tm   l^u(x,y)|    =  h^Ui 2^ li     =  bh   . 

0<r<2h  ^x  I  J 

By    (9.13)    and    (9.20)   x^e   obtain  for    (9.1^!-) 

(9.21)  T^^  TLH   |^u(x,y)<  <  m,+raph   . 

0<r<l     /\_x  ^     '^ 
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It  remains  to  estimate 


'o. 


(9.22)   h2  Z:  Z  l^u(x,y)  |  ^h^TT  eIt-^^^V^)  \ 

^^  Rt<?<R   A^ 

We  recall  that  in  the  above  region 

|u(x+l,y)-u(x,y)  |    <  A:|log  /(x+l^+y^.io^  ^^^+f  \+rrJ^     ..    \.^ 

hence 

f^u(x,y)  I  =    |i  f^u(x;y)  |    <  ^Uuog  s/cx^IFTF-log  s^¥^\  + 

Ax  ^  Ax  i 


/\^x  "(^'^'^     r'^       r'^     r'^+2x+lj 


Transforming  to  polar  coordinates,    i.e.,   letting  x  =  rcosa, 
we  obtain 

M  1 


l|-'^(^.^>l^^[i^-^^=-i^^^* 


-^y       --v^ 


'A^ -  n  ^'ji        -^  ij.nr'^     r^+2r  cos  a+1 

^  Ir  1    2  .  it2^M+l  M       1 

^     r     l\.Ti  r        (r-1) 

returning  to  the  original  plane,  i.e.,  letting  r  =  rh  we  obtain 

I A  „/  „x  I  ^  Ir  1  h  .  Mh^  .  Fih^   T  '■ 
I— -u(x,y)|  <  h^2¥  r"^ "?""'' 2^  * 

Hence  for  1  <  r  <  R  we  have 
and 
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h^  Z  Z  l^u(x,y)  I  <  h2  z  z  M*^^  z  z  I  * 

^x  l<i;<R  '^'^  l<r<R  r'^ 

(9.23) 


^k 


+  h2  2: 


^ 


I'lh 


l<r<R  (r-h)^ 
but 
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lim  h2  ^  Z  ^  F  =  ^  \  dejirdr  =   (R-l)a 
h~>0         l<r<R '^^  ^       '^^  I      T  -^ 

^  .2ti     R 

lim  h^  X  ZI  %  =  Sh  1  dQ^^  =  2iihI4  log(R-l)    =  aph 
h->0         l<r<R  r*^  4       ^  r*^ 


o       1 

2Tt     R 

lim  h2  2:  Z  -^^  =  Mh  f  dQ(-£^  =  27tr-1h[log(r.h)--^]? 
h->0         l<r<R  (r-h)*^  I      ^l(r-h)'^  ^"^  "^ 

=  27il^  log  xiK'-^fTIH  -F^h^  --^  ^3^  • 
Combining  all  above  bounds  we   obtain 

(9.2i!-)  ^     HH   l=^u(x,y)|    <  c^+Cph   , 

Ax 

Since  u(x,y)  =  u(y,x),  in  a  similar  manner  as  above  we 
would  obtain  the  following  bound  for  u(x,y) 

Ay 

(9.25)  h^  Z  Z  |^u(x,y)!  <  c^+c^h  . 

Using  (9.2l|.)  and  (9.25)  as  v/ell  as  the  fact  that  |g(P)|  < 
<  c  5  (derived  at  the  end  of  the  previoias  chapter)  we  obtain 

""  A 

from  (7.10)  and  (7,11)  the  following  boionds  in  G.  for  ==-e,(P) 

A  ^     Ax  ^ 

and  ==-e,(P),  respectively 

Ay  A 

(9.26)  |^e.(P)|  <c"5+0(h^) 

Ax  ^ 

(9.27)  |^e-,(P)|  <  c'54-0(h^)  . 

Ay   ^ 
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10,  A  Bound  for  the  Horaopieneous  Case 

In  order  to  finish  all  the  preparations  for  finding  the 
bounds  for  ===-e  and  %*-e  from  (7»3)  we  still  must  find  a  bound 

A^         Ay 

for  the  solution  of  the  difference  equation  Ah®?  ~  ^  ^^^   ^^^ 
difference  quotients.  We  proceed  as  folloiijs.  By  (7 .5)  we  have 

Ah®2  =  °    ^^  ^h 

epi   =  -e,(P)on  the  boundary  • 

Hence,  from  the  above  and  (8.20),  vie   obtain 

(10.1)  |e2lB  <  \e^{^)  I  <  cJ;'5|log  h|+0(h2)  . 

But  /\u Qg ( P )  -   0  satisfies  the  maximum  principle,  i.e., 
ep(P)  can  take  on  its  maximum  value  only  on  the  boundary,  hence 

(10.2)  le^CP)!  <  c";'5|log  hi 

at  all  interior  points  of  Gj^. 

Next  we  shall  prove  the  following  theorem: 

Theorem.   Let  G  be  a  square  domain,  i.e.,  |x|  <  b;  |y|  <  b;  G' 
--  its  subdomain  and  u(P)  the  set  of  all  lattice  functions  that 
satisfy  the  difference  equations  /wu(P)  =  0  in  G  and  are  uni- 
formly bounded,  i.e.,  |u(P) |  <  A  in  G.   Then  there  exists  a  con- 
stant A»  such  that 

A  A 

|^^-u|  <At    and    |=^-u|'  <A»    in    G»  . 

Ax  Ay 

At  this  point  ^^re  would  like  to  note  that  in  this,  and  only 
this,  chapter  we  shall  denote  the  difference  quotients  of  a 
function  u  by  u^(x,y),  u-(x,y),  u  (x,y),  and  u-(x,y),  i.e., 

(10.3)  u^U,y)   =  Hi2LJil^-JiL2bJlI  .   u.(^-,y)   =  u(x,y)-u(x-h,y) 


5   .    r.t   A        J 


■iZ>y  ^f 


*}:  O.I   'te?-io  ffl 


,i\r 


■■-•■'Ofc 


•  io  X  jJ  ij  Ji,  f> '  t   &  1 1 J    ceo  1 


■iilJ   wo'xl 
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■i'r,r{   -■  i;:-r'T^1 
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(104)   Uy(x,y)  =  u(x,y^-hM(x,y)  .  ^_(^^yj  ^  u(x,y)-u(x,y-h) 
whereas  in  the  previous  chapters  the  above  vjas  denoted  by  —— u 

A  Ax 

and  u,  respectively. 

Ay 

Similarly  we  define  the  second  differences  of  a  function  u 
by: 
(10.5) 


u  -  =  u(x+h.y)+u(x-h.y)-2u(x.77) 


XX 


(10.6)  u  -  =  ^(^»y+^^)+^(^^y-^)-^^(^>y)    . 


yy         h^ 

Hence  the  operator  At-  defined  by  (1.3)  can  be  written  as; 


(10.7) 


Ah^  =  ^xx+^vv   • 


XX  yy 


To  prove  the  above  theorem  we  shall  employ  the  following 
auxiliary  function 

(10.8)   z(P)  =  u^l+C[u2(P)+u2(p^)+u2(P2)+u2(P^)+u2(P^)] 


where 
(10.9) 


I   =  (x2.b2)2(y2.^2j2  . 


C  --  a  positive  constant  to  be  determined  later  and  P, ,  Pp,  P» 
and  Pi  are  explained  by  the  diagram  below. 


P2(x,y+h) 


P^(x.h,y)'  •  -^ 


P(x,y) 


•  P^(x+h,y) 


P|,(x,y-h) 


.^.i.. 


;1.'J1) 


^Ur-Ji, 


lionu' 


-if 


,-;?« 


baa 


+  X>  -' 


CV 
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Ci    ejdi   0  c 


(V*OI) 


>•  ,    ^  V 
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Now  we  shall  show  that  the  function  2(P)  satisfies  the  In- 
equality /\yZ(?)    >   0.  Applying  the  operator  /\,  as  defined  by 
(10.7)  on  both  sides  of  (10,8)  we  obtain: 

(10,10)   Ah^  =  Ah(^x  i)+CAht^^(P)+^^<Pl)+^^^P2^'*'^^(^3^"'^^^V^ 
Let  us  find  first 

(10,11)  Ah^^x  5^  =  ^^x  D  -  •"  (^x  5^ 

XX         yy 

To  compute  the  above  expression  we  make  use  of  the  follow- 
ing formula  which  can  be  easily  verified: 


(10.12)  (fg)    .  =  f g   .  +  f^Sx   "^  ^-S-  +  f     g      . 

XX  XX  X.  X  XX 


Hence 


(10.13)        {v^D       =  u^  J       +   (u^)     I     +   (u^)     i     +   (u^)       5 

^XX  ^XX  XXX  ^xx  ^xx 

Using  (10.12)  once  more  we  obtain  for  (u^)   : 

^  XX 

(10,11;)    (u^)   =  (u^.u  )   =  2u^  u    +  uj^  +  u^   . 

^  XX       ^   ^  XX       ^   XXX     ^^     XX 

2         2 

To  find  (u~:)   and  (u^)  we  make  use  of  the  following  two 
XX       X  - 

formulas,  which  can  also  be  easily  verified: 


•  r       !■•  r 


r-::.;  ;  m^    ev 


C(^ir'^' 
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(10.15)  (fg)x  =  f(Pi)  gx  +  ^x  s 


(10.16)  (fg)_  =  f(P  )  g_  +  f^  g 


X       ^    X     X 


Then 


(10.17)       (u2)^  =  (u^  .  u^)^  =  u^  u^^  +  u^^u^  (P^) 


and 


(10.18)       (u^)   =  (u  .  u  )   =  u  u   +  u   u  (P.) 

^X      ^     ^X     ^XX     XX   ^^ 


Using  (10. ll;),  (10.17)  and  (10.18)  we  obtain  from  (10.1^3) 


=  u2 

XX     ^   XX 


(4   i>  „^  =  ^x  i  ,  ^  f^x  ^xx  -^  ^xx  ^x  (^^1^  J  ix  ^ 


(10.19)  +  tu^  u   +  u   u  (P-,)]  J  + 


XX     XX 


+  [2u  u    +  vT^  +  vT   ]   "^ 

^       XXX     ^^     XX   ~ 


In  a  similar  manner  using  analogous  formulas  for  the  diffe- 
rences  in  y  we  would  obtain  for  (uf  J)   the  following  expression: 

^^'  ^^,^  =  ^K  ^y.  ^  f^x  -xy  -  -xy  -X  ^h^^ly^ 

(10.20)  +  [u^  u  _  +  u  _  u^  (P.  )]  J  + 

xy    xy      ^    y 


+  [2u_  u    +  u^^  +  u*^  ]  J  . 

"^  xyy    "^y    xy- 
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59 
Combining  (10.19)  and  (10,20)  we  obtain  from  (10.1)  the  fol- 
lowing  expression  for  /\,  (u  J): 

A^(u  lb)   =  u  [T   +  J  ]  +  [u  u   +  u   u  (Pn)]  T  + 

'-^h^  X  ^'  X^\'         ^  r       X   XX     XX   X  ^^l'-"  ±x 

XX    yy 

+  [u_  u   +  u   u  (P-)]  £  + 
^   XX     XX   ^    -^     X 

(10.21)  +  [u^  u^y  +  u^y  u^  (Pg)]  iy  + 

+  tu^  u   +  u   u   (P,  )]  I  + 
^  xy    xy  ^   '^    y 

+  [2u  (u    +  u   )+u^^+u^  +^v"*"^^  ^i 
^  XXX    xyy   ^^  xx  ^^  xy 

Remembering  that  one  of  our  assumptions  has  been  that 
^  u  =  0,  i.e., 

(10.22)  u   +  u  _  =  0 

XX    yy 

and  taking  the  difference  quotient  with  respect  to  x  of  the  above 
expression  we  obtain 

(10.23)  u  ^  +  u  ^  =0 

XXX    yyx 

It  is  easily  seen  that 

u  _^  ~  ^  .  ^^^  ^  _  ~  ^ 

XXX    XXX        yyx    xyy 

therefore 

u    +  u    =0 
XXX    xyy 

and  (10.21)  becomes: 

Ai,  [u2  J]  =  u^  [|^  -H  l^J    +  [u^  u^^  ^   u^^  u^  (P^)]  i^  + 


+  [u^  u   +  u   u^  (P,)]  I  + 

^   XX     XX   ^  ^    ~x 

(10.2U)                  +  [u^u^y  -^"xy  ^x  (^2^^  ^y  "^ 

+  [u_  u   +  u   u^  (Pi,)]  i  + 

^  xy    xy  ^  ^    y 


+  (uj  +  u*^^  +  uj  +  u'^^)  5 
^^    XX    ^y    xy 
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To  complete  the  computation  of  /\i,z  we  still  have  to  find 
Ah^^»  Ahu2(P^),  At/(p^)^  Ah^^(P3)  ^nd  Ai/i\)' 

(10.25)  Av,u^  =  (u^)  -+(u^)  -  =  (u»u)  .+  (u.u)  -  . 

— "■  XX    yy      XX     yy 

Applying  (10.12)  to  (10.25)  we  obtain 

(10.26)  Av.^^  =  2u(u  .+u  .)+u^+u^+u^+u^  . 

^        XX  yy   ^  X  -y  y 

Considering  (10.22),  (10.26)  becomes 

(10.27)  Ah^^  =  u^+ut+u^+u§  . 

"■         ^    X    y    ^ 

In  a  similar  manner  as  above  v;e  shall  get  the  follovjing 

o  o  o 

analogous  expressions  for  A^^  (^i)*  Ah"  ^^2^ '  ^h"  ^^3^  ^"^ 

Ah"^(\^»  respectively, 

(10.28a)   Aiy(Pi)  =  "x(Pi)+"|(Pi)+Uy(Pl)+^!(Pl) 
(10.28b)   Ah"^(P2^  ^  u^(P2)+u?(P2)+u^(P2)+uS(P2)  • 
(10.28c)   Ah"^(P3)  =  ul{P^)+ul{?^)+n^^{F^)+nli?^) 


(10.28d)   Ah^^^V  ==  ^x(\)^"^(\)^"y(\^^";(V 


y 

X  '1-   0-   '+   y 


To  obtain  the  upper  bounds  for  some  terms  in  (10,2l|)  we 

have  to  examine  the  function 

i  =  (x2.b2)2(y2.i,2j2  ^ 

Obviously  ^  is  a  continuous  function  for  all  x,y  in  the  closed 
bounded  domain  G:  |x|  <b,  |y|  <b  and  is  tuice  dlfferentiable. 
Let  us  differ  ntiate  J  with  respect  to  x 

(10.30)     ^  =  I^(x2.b2)(y2-b2)2  =  ipc(y2-b2)v^  . 

2      2 
Prom  the  above  vie    see   that  Ipc(y  -b   )    is   continuous   in  the 
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closed  bounded  domain  G;  hence  it  assiunes  its  greatest  value  at 

least  once  and  there  is  a  constant  Q,  such  that  in  G,  max  l|.x  • 

2  ?  x,yeG 

•  (y  -b*^)  =  ©1  and  from  (10.30) 

(10.31)  ll^^l  =  |lpc(y^-b2)  I  s/f  <  ^1  v/l  . 

As  we  have  seen  J  is  a  continuous  function  in  the  closed 
bounded  domain  G  and o J/jx  existsj  hence,  applying  the  mean 
value  theorem  to  J,  we  will  be  able  to  obtain  a  bound  for  the 

difference  quotient  of  J 

*(x+h,y)-J(x,y)  ^±{x+<lb.,y)     ■ 

(10.32)  i^=  H =  ^^—  . 

Considering  (10. 31)  we  obtain 

(10.33)  15^1  <  ©i\/|  . 
In  a  similar  manner  x-je   would  obtain 
(10, 3U)  II-I  <  aiVI  . 

Using  reasoning  similar  to  the  above,  only  with  respect  to  y 
we  shall  obtain  the  following  bounds  for  the  difference  quo- 
tients of  5  in  the  y-direction 

(10.35)      l^^i  <  Qp  \/f    and    ||.  |  <  Q^\/% 

-y    ^  -         y    '^ 

2   2 

where  Pp  =  max  l|.y(x  -b  ) . 

Using  (10,32)  we  shall  show  that  ?  .  is  bounded.  By  (10,32) 


(10,36)   ?  „  = 


—XX 

^5(x+v^h,y)  aJ(x+T''h-h,y) 


^Hk-,  .y) 


1'  3x2  ■ 

Differentiating  (10.30)  with  respect  to  x  xie  get 


(10.37)        -^   =  il-(y^-b2)  v'^+8x2(y2-b2)2  . 

9x 

hence  o)  J/3x  is  a  continuous  function  in  the  closed  bounded  do-^ 
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main  G;  |x|  <  b;  |y|  <  b;  therefore  it  achieves  its  greatest 
value  at  least  once,  i.e. 

(10.38)  I p — I  <  max  ^ —  =  Q.  , 

3x       x,yeG   'cx'^^      -^ 

Using  the  above  we  obtain  from  (10. 36) 

(10.39)  IS^-I  <  63  . 

In  an  analogous  manner  we  obtain  the  following  estimate 
for  5  . 

(io.i;o)  li  _|  <  0. 

yy    ^■ 

where  « 

'rl(x,y) 


Qi  =  max  — 5 —  • 

^   x,y8G   By*^ 

Combining  (10,39)  and  (lO.i^O)  we  obtain 

,   V  XX     YY      J)   M-     P 

(lO.i^l)  ^^ 

Choosing  9  such  that  0  =  max(©,  ,©p,e-,9|-)  vie  shall  have  the  fol- 
lovjing  bounds  for  the  various  difference  quotients  of  J, 

Using  the  above  bounds  and  the  fact  that  for  any  two  numbers  a,p 

2  2 
(10.1+3)  ap  <  S^l^  <  a^+p^ 

we  obtain  for  any  e  >  0 

2  2 


••) 
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^  '\l 
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XX^X 
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(10.i|5)  ILV^  -I    <  -f +e2a2|;u2_ 
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x'^l'  ^„2^2t.,2 


(10.1+8)  lix^x^^i^^xxl  < -^^-~-+e''o''i:u; 


e 
2 


(10.U9)  ISx^x(^3^"-=l    <  -^S-^+e^^^iu^- 


XX  e2  XX 


(10.50)  \ly\(fz)\y\    5  ^2-^-e2925u2y 

(10.51)  |Ju^(P.)u  .1    <  ^r-^+e2o2|u2__   , 

'-y  -^     ^     xy     -       e2  ~  xy 

Considering  the  above  bounds  we  conclude   that    (10,2i;)   re- 
duces to 

Ah(^x2)    >   (l-2c^e2)|(u2   +u2.+u2   +u2.)-(iL+o)u2- 

(10.52)  ^^  "^^       ^ 

.^[u2(P^).u2(P2).u2(P3).u2(P^)]    . 

If  we  add  to  both  sides  of  the  above  inequality  the  following 
positive  expression 

and  substitute  on  the  right  side  the  computed  values  for  At-U  « 

A^u^iP-^),   Ah^^(P2^'  ^h''^^^3^  ^'''^   Ah^^(\)  ^^°m   (10.28a)  to 
(10.28d),  respectively,  we  obtain 
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61; 

>  (l-2e2a2)i(u2   +u2.+u2  +u2.)-(C-iL.o)u2+(c--L) 
•   [u2(p^)+u2(P2)+u2(p  )+u^{P.  )]+C[u?+u?(P3_)  + 

2  /  n    \-L,,2/r>     \j.,,2/T,    \  1  in  r  ,,2  ,  .  2  /  T-,    \  ,  ,,2 , 


^:(P2^'^"-^^3^"*"^-(^[|.^^"^^f^y'^"/^l^"*'^/^2^^ 


+ 

X       '■        X       -^        X 


+  u^(p  )+u^(P.  )]+C[u!+u5(P^)+u?(Pp)+u!(PJ  + 
y3yi|  yyly2y3 

According  to  (10.10)  we  Identify  the  left  side  of  the  above 
inequality  as  A.v.^>  ^^^  ^^  "®  select  e  and  G  so  that 
(10.5i|-)      e^e^  <  i    and    C  >  -^+0 

which  also  iraplles  that  C-l/e  >  0;  we  conclude  from  (10,53) 

(10.55)  Ah^  >  0  • 

Writing  out  the  above  operator  explicitly  we  get 

z(P)  <  ^[z(P^)+z(P2)+z(P2)+2(P^^)] 
at  all  lattice  points. 

Consequently  z;(P)  can  have  no  maximum  at  an  interior  point 
of  any  set,  although  it  may  hnve  a  minimum.  Hence  the  maximum 
value  of  z(P)  must  occur  on  the  boundary.   But  by  (10,9)  J  =  0 
on  the  boundary.  Therefore  in  the  whole  square  |x|  <  b;  |y|  <  b 

0  <  2(P)  <  Zg 

where  z^  denotes  the  value  of  z  on  the  boundary.  Considering 
(10,8)  vie   obtain 

0  <  z(P)  <  5CA^  , 
Since  the  second  term  of  (10,8)  is  non-negative,  we  con- 
clude that  for  P  e  G* 
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(10.56)  u2(P)  <  ^ 

1 


where  for  J  we  have  to  take  Its  lovrer  bound. 

Let  us  denote  by  d  the  distance  between  G  and  G';  then  from 

i  =  (x2-b2)2(y2 .1,2)2  ^  (;c-b)2(x+b)2(y-b)2(y+b)2  , 

taking 

|b-x|  >  d 

|b+x|  >  d 
as  well  as 

|b-y|  >  d 


|b+y|  >  d 

nince  d   : 

side  of  (10,56)  for  (A')^  and  thus  the  following  estimate  for 


we  can  majorize  J  by  d  ,  Since  d   >  0  we  can  take  the  right 


u  is  established 
^                             v/5CA 

(10,57)            "x  -1^'  =  ~Tir~  • 

In  a  similar  manner  we  would  obtain 

(10.58)     |u  1  <A«  ,   u.  <A»   and 

u- 

y 

<  A'  . 

Since  Qpi?)    satisfies  all  the  assuiiptions  of  the  theorem 
proved  above,  with  A  =  c "51  log  h|  we  conclude  that  in  the  sub- 
domain  G'  we  shall  have 

v/5ga 

(10.59)  1(^2^x1  :^^'  =  — —  =  °o^'^°S  h\ 

(10.60)  l^^a^y'  "  °o^'^°S  h\ 

(10.61)  |(ep)-|  <  c  5Jiog  h| 

^  X     " 

(10.62)  1(62)-!  <  CoSilog  h|  . 
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11,  Estimates  for  (A/A^)^  ^^^  (A/A^)® 

Before  using  the  subsidiary  estimates  of  the  previous  chap- 
ters  to  obtain  the  bound  for  ^==~e  and  ^=^e,  we  recall  that 

A^         Ay 

(11.1)  ^-e  =  gl^-L^^^Q.l;?  -^,  >XL 

(11.1)  ^^e  ^^ 

(11.2)  ^e  =  e(x,y+h)-e(x,y-h) 

(11.2)  ^^e  ,^ 

and 

(11.3)  ©  =  ®i'^®2  • 

Therefore  talcing  difference  quotients  of  both  sides,  we  get 

A       A    ,A 

Ax         Ax       A^ 
hence 


(11.^)  I^el    <   |4-eJ  +  |4 


and  similarly 
(11.5) 


Ax  ■  "  'A^  "    Ax'^ 


Ay        Ay       Ay  • 


A  A 

However,  in  the  estimates  for  ===-6^  and  ==^e«  as  well  as 

A  A  A^  A^ 

for  =^e,    and  ^=^6^,  given  by   (9.26)7  (9.2?) ,    (10.59)    and 

Ay  ^    Ay  ^ 

(10,60)  the  forward  differences  have  been  taken,  while,  accord- 
ing to  (11.1)  and  (11,2)  centered  differences  are  required. 
Therefore  having  at  our  disposal  the  estimates  for  baclward  dif- 
ferences as  \-jell,  and  noting  that 

^e  =  |(e  -i-e.) 
Ax    2  X  X 
and  similarly 


.4™e  =  i(e  +ej  . 

Ay      2    ^    y 


Finally  we  t^jould  obtain 


(^:iv-...:?  ^--^^ 


*i,Ux     ViJ,J<^'\i 
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Ax 


(11.6)  |^e|  <  C^5|log  h|+0(h^) 


and 


Ay 


(11.7)  l^e|  <  C^5|los  h|+0(h'^)  . 


12 .  R e fie c tion  Principle 

In  previous  chapters  we  found  the  estimates  for  e,  e  and 

A  ^"^ 

==-Q  where  e  satisfied  the  following  equation 

Ay 

Ah®  "^   S(^^   ^^  ^^®  interior  of  G^^:  |x|  <  1,  |y|  <  1 
e  =  0     on  the  boundary  . 

A  A 

Hovjever.  our  estimates  for  ==--e  and  ==-e  do  not  hold  suffl- 

A^      Ay 
ciently  close  to  the  boundai-'y.   To  make  them  hold  up  to  the 

boundary  we  shall  proceed  as  follows. 

Suppose  we  are  given  a  function  e(x,y)  vjhich  is  the  func- 
tion e  described  above,  and  lot  us  define  a  function  E(x,y)  in 
the  domain  D;  |x|  <  3,  lyl  <  3  such  that 

A^E  =  G(P)   in  D 

E  =  0     on  the  boundary  of  D 
and 

(12.1)  E(x,y)  =  e(::,y)  in  the  interior  of  rectangle  I: 

-l^x<l  ;  -1<7<1 

(12.2)  E(x"''",y'''")  =  E{-2-x,y)  =  -e(x,y)  in  the  interior  of  a 

rectangle  I„:  -3  <  x""'  <  -1  ;  1  <  y"'  <  3 

(12.3)  E(x'^,y'"')  =  E(x,2-y)  =  -e(x,y)  in  the  interior  of 

rectangle  I^:  -1  <  x'"  <  1  ;  1  <  y"  <  3 
(I2.I4-)  E(x'"',y'"")  =  E(2-x,y)  =  -e(x,y)  in  the  interior  of 

rectangle  I„:  1  <  x""  <  3  ;  -1  <  y"  <  1 


bnsi 


m\ 


i'  ^'^'^ -f "  (■ ' ■'•■■-' 


.....    ,'.;f.    is  -'-r.    ? 


;,  i.rJf'Y     <^r*:i'     •■'•  J 


t  '-        r 

1    I    •  - 


.j.      .-,.r.+      p. 


'^■.C.f.^l-^M,---,  ■ 


.    .J 


XIu*?- 


6f] 

(12.5)  E(x",y"'^)  =  E(x,-2-y)  =  -e(x,5)  in  the  interior 

rectangle  I^:  -1  <  x"''  <  1  ;  -3  <  y*''  <  -1  • 

But  in  I 
a 

(12.6)  AhE(x»y)  =  -Ah®(P)  =  -s(P)  =  G(P) 

and  similarly  in  I^,  I^  anci  I^  we  see  that  G(P)  =  -g(P);  it  is 
obvious  that  G(P)  =  g(P)  in  I  (see  Fig.  2). 

Let  us  take  noi\r  a  point  B  on  the  boundary  of  I  and  I  and 
compute  the  value  of  E  at  B  (denoted  by  E„)  by  considering  the 
points  numbered  la,  2a,  3a,  l;a,  and  B  in  I.  and  points  numbered 
1,  2,  3,  U-   and  B  in  I  (see  Pig,  3)  and  writing  out  /VE  =  G(P) 
explicitly  for  those  points  vre  obtain: 

(12.7)  ^l^'-^Za^^3^'-^^-^k^   =  V 

(12.8)  Ej^+Eg+E^+Eg-lj-E^  =  G^  , 

But  by  definition: 


(12.9) 


as  well  as 


(12.10) 


%  =  ^1 

^2  =  ^2 

E3  =  e3 

\  =  e^^     and 

%- 

=  % 

%a  =  -^1 

E2a  =  -^2 

^3a  =  -^3 


\a   =  '%     ^^<3  \a   =  -%  ' 
therefore  (12.7)  and  (12.8)  become 

.e^-e2-e3+Eg+l|e^  =  g^ 

e^+e2+e3+E3-l^e^  =  -g^  ; 

adding  the  above  equalities  v;e  get 


?.h 


r.   ri    r--::f' 


"  r  V 
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Sj^ 


■'^!   r^i^'*^ 


Tr-.^.f -pf  ^ .-     x.-j-.j. 
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•?^-*'*^,"V      ^^"4  +1 ' 


li  ,• 
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2E  =  0 

(12.11)  ^ 

Eg  =  E(-l,y)  =  0  . 

In  a  similar  manner  talcing  corresponding  points  in  I  and 

Ivi  I  and  I^,  I  and  I,  we  irould  obtain 
D'        c '        d 

(12.12)  E(-l,y)  =  E(+l,y)  =  E(x,l)  =  E(x,-1)  =  0  . 

Now  we  are  able  to  compute  the  value  of  G  on  the  boundary 
of  |x|  <  1,  |y|  <  1.   Using  points  numbered  l^   and  i|a  in  I  and 
I  respectively  and  three  boundary  points  and  applying  Ai-E  = 
=  G(P)  to  them  vie  obtain 

(12.13)  Ej^^+E^+Eg+Eg-UEg  =  G3  , 
but  by  definition 

and  from  (12.11) 


£3  =  0; 


hence    (12.13)   becomes 


■•%^%  =  % 


^12. 111.)  Q^  =  G(-l,y)    =  0   . 

Similarly  taking  points  in  I  and  I,  ,  I  and  I  ,  I  and  I,  we 
would  obtain 
(12.15)    G(-l,y)  =  G(l,y)  =  G(x,+1)  =  G(x,-1)  =  0  . 

It  is  clear  that  what  we  did  above  vras  obtained  by  reflect- 
ing the  rectangle  I  about  lines  y  =  +1  and  x  =  +1.  Repeating 
this  process  once  more,  i.e.,  reflecting  I  about  y  =  +1,  and 
I  about  y  =  il,  0^  I^  about  x  =  +1  and  I,  about  x  =  +1  in  a 
similar  manner  as  before  we  ijould  obtain  the  following  rectangles 
which  will  be  denoted  by  II_  and  II""  and  II.  and  II-;!-  or  II,  and 
11^  and  11^  and  II'^",  respectively  (see  Pig.  2). 
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Using  the  same  methods  as  above  it  Is  easily  verified  that 
in  the  interior  of  11^;  -3  <  ::'"""'  <  -1;  1  <  y''"""  <  3,  E(x""''"  v""'")  = 
=  -E(x'^,y")  =  e(x,y),  and  G(P''""')  =  -G(P''")  =  g(P)  and  that 
E(x''^"",l)  =  E(-l,y"""')  =  0  as  x.-ell  as  G ( x"""''" ,  1 )  =  G(-l,y""''"')  =  0; 
in  the  interior  of  11^":  -3  <  x'"""'  <  -1,  -3  <  y"""'^  <  -1,  E(x""^",y'''""') 
=  -E(x'''',y'')  =  e(x,y)  and  G(P'""''')  =  -G(P''')  =  g(P)  and  that 
E(x'""",-1)  =  E(-l,y'"""')  =  0  as  well  as  G(x"'^,-1)  =  G(-l,y'"""')  =  0; 
in  the  interior  of  I,  :  1  <  x''""  <  3;  1  5  y'"""'  <  3,  E(x'"""',y '""■')  = 
-E(x''^y'"')  =  e(x,y)  and  G(P'''""'')  =  -G(P'''")  =  g(P)  and  that  E(x"^^,l) 
=  ECljy"^^')  =  0  as  well  as  G(x''"%l)  =  G(l,y'"")  =  0;  In  the  In- 
terior of  11^":  1  <  x'"""'"  <  3,  -3  <  y'"""'  <  -l,  ECx'^'^y''^")  = 
-E(x'"'',y*)  =  e(x,y)  and  G(P'''~'')  =  -G(P")  =  g(P)  and  that 
E(x^'"'',-1)  =  E(l,y'"""')  =  0  as  v/ell  as  G(x'''"",-1)  =  G(l,y''^"')  =  0. 

Therefore  we  conclude  that  the  function  E(x,y)  defined  in 

the  domain  D:  |x|  <  3,  |y|  <  3  satisfies  /Ve   =  G(P)  in  D  and 

E  =  0  on  the  boundary  of  D  and  that  E  as  vjell  as  G  is  zero  on 

the  lines  x  =  +1,  y  =  +1  and  othervjise  E(P"")  is  either  +e(P)  or 

-e(P),  as  well  as  G(P'"')  is  either  +g(P)  or  -g(P).  Hence  in  the 

Interior  of  D 

|E(P'")1  =  |e(P)|  <  6 
and 

|G(P-')|  =  |g(P)|  <  0^5  . 

Therefore  i/e  are  able  to  apply  all  the  theory  developed  in 

chapters  7-11  to  the  function  E  and  conclude  that  in  the  subdo- 

main  D';  |xl  <  2,  |y|  <  2,  the  follovjing  estimates  hold: 

A 

(12.16)  |-^i-E|  <  C^5|log  h| 
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(12.17)  |~-El    <  C   5|log  h|    . 

Ay         ° 

But  E(x,y)    =  e(x,y)    In  Q,  :    |x|    <  1,    |y|    <  1;   hence    the 

A  A 

following  estimates  for  — — e  and  ==^-6  hold  up  to  the  boundary 

A^         Ay 

|^e|  <  G  5|los  h| 

l=f"6|  <  C  Sjlog  h|  , 

Ay 


4; ;  ■*.':•'.  i 
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